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Inelastic Column Theory 


F. R. SHANLEY* 
Lockheed Aircraft Corporation 


SUMMARY 


The action of a column in the plastic range is analyzed on the 
basis that bending may prageed simultaneously with increasing 
axial load. This leads to a new column formula that includes 
both the tangent-modulus (Engesser) and the reduced-modulus 
(von Kérm&n) formulas. It is shown that bending begins at 
the tangent-modulus load and that the column load increases 
with increasing lateral deflection, approaching the reduced- 
modulus load as a limit if the tangent modulus is assumed to 
remain constant. 


INTRODUCTION 


I A RECENTLY PUBLISHED PAPER,! the author 
stated that the reduced-modulus (or double- 
modulus) theory is not correct for predicting the 
maximum load up to which a perfect column will 
remain straight. This is because it is possible for the 
column to bend simultaneously with increasing axial 
load. Under such conditions it is possible to have 
bending without introducing any strain reversal, 
upon which the reduced-modulus theory depends. 
In reference 1 it was stated that the column will begin 
to bend as soon as the axial load exceeds the value 
predicted by the tangent-modulus (Engesser) theory. 
It appeared likely that the Engesser load could be 
exceeded but that the reduced-modulus load could not 
be reached. In this paper it will be shown that for an 
idealized simplified column, this is actually true. 

The three basic column formulas may be written as 
follows (assuming pin ends and zero eccentricity): 


(Euler) P, = r°EI/L? (1) 
(Engesser) P, = r°E,I/L? (2) 


(reduced P, = (3) 
modulus) 


where 


Presented at the Structures Session, Fifteenth Annual Meeting, 


1A.S., New York, January 28-30, 1947. 
* Division Engineer, Engineering Research. 


= critical load 

= moment of inertia of column cross section 

column length 

= Young’s modulus (slope of stress-strain dia- 
gram in elastic range) 

E, = tangent modulus (local slope of stress-strain 

diagram in inelastic range) 
E, = reduced modulus (an effective value lower 
than F and higher than £,) 


Eqs. (1), (2), and (3) differ only in the value used 
for the effective modulus of elasticity. Since the 
Euler equation applies only in the elastic range, the 
problem of column action in the inelastic (plastic) 
range centers around Eqs. (2) and (3). An excellent 
summary of the history of these two theories is given 
in reference 2, from which the following is quoted: 

‘“‘What is here called the double-modulus theory has 
frequently been called also the Considére-Engesser 
theory and K4rm4n’s theory. Many competent engi- 
neers are mistaken as to the origin of the theory, and 
a brief account of its development will not be out of 
place. In 1891 there was published a memoir included 
as an appendix (annexe) to the proceedings of the 
Congrés International des Procédés de Constructions, 
held in Paris from the 9th to the 14th of September, 
1889, in which A. Considére pointed out that, as an 
ideal column stressed beyond the proportional limit 
begins to bend, the stress on the concave side increases 
according to the law of the compressive stress-strain 
diagram, and the stress on the convex side decreases 
according to Hooke’s law, and that therefore the 
strength would be given by 

P = 
in which £ is a modulus, the value of which lies between 
the modulus of elasticity and the tangent modulus. 
Considére realized that E was a function of P/A, the 
average stress in the column, but did not go further 
than to point out this fact. 
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a> “a V—EULER EQUATION 
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Fic. 1. Column theories and test data (Alcoa). 
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40 60 80 100 
SLENDERNESS RATIO, 4 


Column theories and test data (Van den Broek). 


Fic. 3. Pin-end column test arrangement. 


“Earlier in 1889 Fr. Engesser presented his tangent- 
modulus theory, and in 1895 Félix Jasinski (“Noch ein 
Wort zu den Knickfragen,’ Schweizerische Bauzeitung, 
Vol. XXV, No. 25, p. 172, June 22, 1895) pointed 
out that this theory was not correct and called attention 
to Considére’s work. He stated that at that time it 
was impossible, however, to determine theoretically 
the form of the function £. Thereupon Engesser 
(‘Ueber Knickfragen,’ Schweizerische Bauzeitung, Vol. 
XXVI, No. 4, p. 24, July 27, 1895) acknowledged the 
error in his original theory and replied that the possi- 
bility of determining E theoretically was in no wise 


out of the question, and he determined it in the genera] 
Nothing further was done apparently 
until Karman (‘Untersuchungen iiber Knickfestigkeit, 
Forschungsarbeiten, Nr. 81, 1910) presented the theory 
again, added the actual evaluation of £ for the rectangy. 
lar cross section and the idealized H-section (consisting 
of infinitely thin flanges and negligible web), and gave 
the theory new life by making a series of careful tests 
designed to afford a check on the theory. Since then 
E has been evaluated for other cross sections by a num- 
ber of writers.” 

Although the reduced-modulus (or double-modulus) 
theory has long been accepted as the exact theory of 
column action, the simpler tangent-modulus theory 
has been found to be much easier to use, since it js 
not affected by the shape of the cross section. Since 
it gives lower critical loads than the reduced-modulus 
theory, it is also preferred by engineers on the basis 
of safety. Finally, test data indicate that the actual 
buckling loads are usually closer to the tangent-modulus 
values than to the reduced-modulus values. A typical 
example is shown in Fig. 1, which is reproduced from 
reference 3 (tests by the Aluminum Research Labora- 
tories). Reference 4 also shows test data that agree 
with the Engesser theory. Fig. 2 is based on Fig. 4 
of this reference, which includes a number of tests 
made by Professor Van den Broek. 


Test DATA 


In an effort to obtain a clue, the author had a simple 
test performed by the Lockheed Research Laboratory. 
A pin-ended 24ST aluminum-alloy column of rectangu- 
lar cross section (2 by 11/4 in.) was designed to fail in 
the plastic range and was equipped with electric strain 
gages on opposite sides, at the mid-section (Fig. 3). 
(Details of the test are omitted here but are reported 
in reference 5.) 

The results obtained are plotted in two different 
ways in Figs. 4 and 5. 

Fig. 4 shows that the strain distribution remained 
substantially constant up to a load of about 40,000 lbs. 
and then gradually shifted, indicating bending. Up 
to a load of about 87,000 lbs. there was no strain re- 
versal. 
place. At the maximum load of 92,500 Ibs., con- 
siderable strain reversal had taken place and the axis 
of rotation of the cross section appeared to be at about 
one-third of the width of the column. 

Fig. 5 shows the strain on opposite sides of the 
column plotted against load. The slight divergence 
corresponds to an initial eccentricity of about 0.001 in. 
If the column were perfect and remained straight up 
to the reduced-modulus load, the two curves would 
coincide and would be represented by the dashed 
line OA. (The upward curvature reflects the shape 
of the compression stress-strain diagram.) Point B 
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Fic. 4. Strain distribution as determined in a column test. 
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Fic. 5. Strain on opposite faces of column from test data. 
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represents the upper limit as predicted by the tangent- 
modulus theory. 

Fig. 5 shows that if the column were to remain 
straight up to the reduced-modulus load there could 
be no strain reversal below that load. What, then, 
can supply the extra effective value of E needed to 
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Fic. 6. Alternative types of strain distribution across column 
cross section. 
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Fic. 7. Type of strain distribution needed to permit loading 
beyond tangent-modulus load. 


prevent buckling beyond the tangent-modulus load? 
The obvious answer is that the column cannot remain 
straight beyond the tangent-modulus load; there 
must be a definite amount of strain reversal as soon 
as the load is further increased. This should cause 
the two curves to separate at point B, one starting 
downward and the other upward. 

It can-now be seen that in the derivation of the 
reduced-modulus theory a questionable assumption 
was made. It was assumed, by implication at least, 
that the column remains straight while the axial load 
is increased to the predicted critical value, after which 
the column bends, or tries to bend. Actually, the 
column is free to “try to bend” at any time. There 
is nothing to prevent it from bending simultaneously 
with increasing axial load. Under such a condition it 
is possible to obtain a nonuniform strain distribution 
without any strain reversal taking place. The differ- 
ence between the two assumptions is shown diagram- 
matically in Fig. 6. 

Fig. 6b, however, still represents a paradox. There 
is no strain reversal indicated; hence, the value of 
E, must apply over the entire cross section; therefore, 
the column load cannot exceed the tangent-modulus 
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Fic. 8. Alternative types of strain progression. 


value. If the load is to go any higher, some strain 
reversal must begin at the tangent-modulus load. 
The picture might then look something like Fig. 7, 
in which each succeeding increment of loading beyond 
the tangent-modulus load causes some additional strain 
reversal. The fact that this picture resembles the 
actual distribution shown in Fig. 4 is significant. 
(NoTE: These statements do not necessarily apply 
to the entire /ength of the column, since the ratio of 
bending moment to axial load varies from end to 
end.) 

Another way to visualize the two different theories 
is shown in Fig. 8, in which the points represent condi- 
tions on opposite edges of the column. In the reduced- 
modulus theory it is assumed that the points stay to- 
gether until the critical load is reached (Fig. 8a). Fig. 
8b shows that they may separate without involving 
any strain reversal. Since the average stress in Fig. 
8b is greater than that corresponding to the tangent- 
modulus theory, it represents an impossible condi- 
tion. Fig. 8c shows what would have to happen if 
the average stress were to exceed the tangent-modulus 
value. 

On the basis of the foregoing reasoning the author 
predicted, in reference 1, that (a) bending will begin 
as soon as the tangent-modulus load is exceeded; 
(b) the maximum column load will be reached some- 
where between the loads predicted by the two theories. 


MATHEMATICAL ANALYSIS 


In order to prove the last statement, for at least 
one case, the problem has been greatly simplified by 
adopting a suggestion made by E. I. Ryder, of the 
Civil Aeronautics Authority. This consists in working 
with a two-legged hinged column in which the hinge 
consists of a unit “‘cell’’ formed from two small axial 
elements. As shown in Fig. 9, the two legs of the 
column are assumed to be infinitely rigid. If the 


dimensions of the unit cell are sufficiently small with 
respect to the column length, L, it can be assumed 
that there is a simple hinge action about the center of 
the cell. This device reduces the problem to ele. 
mentary form by eliminating the mathematical work 
involved in integrating over the cross section and over 
the length of the column. 

As shown in Fig. 9, the two elements of the column 
cell are assumed to have deflected in opposite directions 
through the distances e,; and é,, which may be regarded 
as the strains that occur after the column starts to 
bend. If e; and e, are assumed to be equal and opposite 
in sense (as shown in Fig. 9), pure bending is indicated, 
This restriction will not be applied, however; e; and ¢ 
may have different values, indicating combined bend- 
ing and variation in axial load. 

The critical load for such a column is easily deter- 
mined by equating external and internal bending 
moments. The lateral deflection, d, is first expressed 
in terms of strain: 


a- +8) @ 


The external bending moment at the hinge is 
M, = Pd = (PL/4)(e: + e2) (5) 
The axial load in each flange element, due to bending, is 
P, = @F,(A/2) 


Note that E; and £, indicate the value of E which is 
considered to be effective for each flange element. 

The internal bending moment (about the hinge 
point) may be expressed as 


[ 
INFINITELY 


Fic. 9. Simplified two-flange column. 
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(ek: + ak) (7) 


Equating internal and external bending moments 
(Eqs. (5) and (7)], 
(PL/4)(é1 + = + 
from which 
P = (A/L)[(e:E, + e2E2)/(e1 + &)] (8) 


Now in either the Euler or Engesser equations E; = 
FE; Eq. (8) therefore reduces to the following two 
equations: 


(Euler) P, = AE/L (9) 


(Engesser) P, = AE,/L (10) ° 


If it is now assumed that element (1) undergoes 
increasing compressive stress while element (2) has a 
decreasing compressive stress, EZ, and FE, may be 
replaced by E, and E£, respectively. 

Letk = E/E,. Then 


= E, 
E, = kE, 
Substituting in Eq. (8), 
P = (AE,/L)[(¢: + + (11) 


This is the same as the Engesser formula except 
for the added term, which will now be further exam- 
ined. 

From{Eq. (4) 

+ = 4d/L 
and 
a = 4d/ L- (12) 
Substituting these values in Eq. (11), 
AE, L [(4d/L) — + ke2)] 


4d 
_ AE, 
pa E | (13) 


Another expression for P may be obtained by assum- 
ing that, after the tangent-modulus load is reached, 
the column load continues to increase. This increase 
is given by the difference between the element 
loads P, and P2, which can be expressed, from Eq. (6), 
as 


AP = P, — = @E,(A/2) — /2) 
Substituting for and 
AP = (A/2)E,(e, — kes) 
Substituting for ¢:, 


A\ _ (4d 


This value should be added to the tangent-modulus 
load to obtain the total value for P. 


Ap | 4d _ 
P = P,+ AP L 


From Eqs. (13) and (15) the following equation can 
be set up: ; 


Ne = (16) 


From this, 
e, = 8d?/L[k — 1 + 2d(1 + k)] (17) 
Substituting this value of e, in Eq. (13), 


L 8) 


This may be reduced to 


Gl) 


If the ratio E,/E is represented by the usual symbol 
7, Eq. (19) becomes 


R = P/P, (20) 
Then, from Eq. (10), 


In Fig. 10, R is plotted against lateral deflection for 
two different values of r. [NoTe: Eqs. (18) to (21a) 
do not apply when £, = 0, since the limiting column 
load is then determined by the stress at which this 
occurs. ] 
It should be noted that Fig. 10 corresponds to the F 
usual plot of column load against lateral deflection. If 
the tangent modulus had been assumed to decrease with 
increasing stress (as it usually does), the curves would 
rise to some maximum value and then start downward 
again. It is interesting to note that in the very short 
column range, where the tangent modulus approaches 
a constant value, test points often lie closer to the re- 
duced-modulus curve than to the tangent-modulus 
curve. 
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Fic. 10. Variation of column load with lateral deflection, 
assuming E; constant. : 


COMPARISON WITH REDUCED-MopULUS THEORY 


The equation for the critical column load will now 
be derived on the basis of the assumptions originally 
used in the reduced-modulus theory. It will be as- 
sumed that the column remains straight up to the 
critical load P,, after which it bends. The derivation 
proceeds as before, up to and including Eq. (13). Now, 
instead of assuming that there can be an increase in 
load, AP, it will be assumed that AP = 0. Then 
P, = P, and Ee, = Ee, from which 


@ = &(E,/E) = e/k (22) 
Substituting for e, from Eq. (12), 
= (1/k)[(4d/L) — 
= (4d/L)[1/(1 + k)] (23) 
Substitution of this value in Eq. (13) yields 
P, = (AE/L){1 + (24) 


Eq. (24) obviously represents the limiting value of 
Eq. (19) as d approaches infinity. This proves that, 
for the simplified column, the reduced-modulus theory 
gives the limiting value for the column load as the 
lateral deflection approaches infinity, assuming that 
the value of the tangent modulus, E,, remains constant. 

The actual upper limit for the column load will de- 
pend largely on the manner in which E, varies with 
increasing strain, as shown by von K4rman.® It is 
important to note, however, that, even for a perfect 
column, lateral deflection must take place as soon as 
the tangent-modulus load is exceeded and that the 
load predicted by the reduced-modulus theory can never 
be reached, even if there is no dropping off in the 
tangent modulus with increasing strain. ~ 

Eq. (19) represents, for the simplified column and 
for a constant value of E,, the complete theory of 
column action. It includes the Euler equation (when 
k = 1) and the Engesser (tangent-modulus) equation 
(when d = 0) and approaches the reduced-modulus 
equation as a limit (asd—~ ©). It is suggested that 


an equation of this type be called the column equilibrium 


equation to avoid any question as to the definition of 
buckling or instability. 


VARIATION OF STRAIN WITH LOAD 


The equation already derived may be used to find 
equations for the strains in the two column elements, 
From Eqs. (21), (12), and (17) the following equations 
may be derived: 


k—R (25) 
2 R-1 q 


From these equations it should be possible to obtain 
a graphical picture of the variation of strain against 
column load, to compare with the experimental curve 
of Fig. 5. Since Eqs. (25) and (26) apply only after 
the column has started to bend, at the load P,, it is 
necessary to determine the average compressive strain 
at the load P;. The stress, obtained from Eq. (10), is 


P,/A = E,/L 


The actual strain corresponding to this stress would 
be obtained from the stress-strain diagram of the 
material. If the stress were in the elastic range, the 
strain would be 


e, = P/AE = E,/EL = 1/kL (27) 


This would be the value of the strain at R = 1. 
(For most materials the strain would actually be con- 
siderably higher, but this does not affect the general 
results.) The additional strain in each element, 
beyond R = 1, will be put in terms of the strain ¢, 
giving 


4m (29) 
k—1 


These values are plotted in Fig. 11, which represents 
a typical variation of strain with column load, for a 
value of k = 1.333. 

Although Fig. 11 is based on extreme simplification 
of the problem and does not give a true picture of typical 
conditions, it has a general resemblance to Fig. 5, 
which was obtained from actual test data. It is 
interesting to note that on the concave side of the 
column the compressive strain increases rapidly after 
the tangent-modulus load is reached, while on the 
convex side the strain starts to decrease rather slowly. 
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The rapid increase of compressive strain which is 
required to obtain additional load beyond P, will 
usually involve a large reduction in E,. This indicates 
that in most actual cases the maximum column load 
will exceed the tangent-modulus load by only a slight 


amount. 


CONCLUSIONS 


Although the foregoing analysis is based on a hypo- 
thetical column that bears little resemblance to an 
actual column, extension of the theory to the more 
general case is largely a matter of mathematics. The 
fundamental principles of column action in the in- 
elastic range will not be changed by such generaliza- 
tion. The following conclusions can therefore be 
drawn from the’ analysis of the simplified column: 

(1) The tangent-modulus (Engesser) formula gives 
the maximum load at which an initially straight, 
centrally loaded column will remain straight. ° 

(2) The column load may exceed the tangent- 
modulus load but cannot be greater than the reduced- 
modulus load. (The latter statement has not been 
proved for the general case.) 

(3) Loading beyond the tangent-modulus load will 
cause bowing, which will produce permanent bending 
deformation (eccentricity). 

(4) There will be some portion of the column cross 
section for which the stress will never exceed the 
tangent-modulus stress. 

(5) After the tangent-modulus load is exceeded, 
the compressive strain over a portion of the cross 
section will increase much more rapidly than the average 
strain. 

(6) For most engineering materials the decrease 
in tangent modulus with increasing strain will limit the 
amount by which the column load may exceed the 
tangent-modulus value. 

(7) The idealized simplified column treated in this 
paper cannot sustain the reduced-modulus load unless 
the column is laterally supported while the axial load 
is being applied. 


o 
3.0 3 31 
25 
[ie 
20 
e zi flo 
t jw 
« 
e 
° 2 4 & & 19 12 
5 
t 


Fic. 11. Calculated strain variation with column load for a 
hypothetical column. 


(8) When supported against bowing during axial 
loading, a column in the inelastic range will sustain a 
greater axial load (with support removed) than when 
loaded without such support. 

(9) The tangent-modulus (Engesser) equation 
should be used as a basis for determining the buckling 
strength of members in the inelastic range. 

In conclusion, it would seem fitting to repeat a 
statement made by von Karman® to the effect that 
“Engesser appears to have been the first to recognize 
properly the nature of ‘inelastic buckling.’’’ It is 
interesting that aircraft engineers, in seeking greater 
accuracy in the inelastic range, have gone back to the 
formula that was first suggested by Engesser over 50 
years ago. 
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Discussion 


Dr. Theodore von Karman: ‘The theory of buckling of a 
column subjected to axial load beyond the limit of elasticity 
has been discussed by many authors in the last three decades. 
Most of the objections, however, were not valid. Now Mr. 
Shanley comes along with an objection that is worthy of con- 
sideration. 

“Both Engesser’s and my own analyses of the problem were 
based on the assumption that the equilibrium of the straight 
column becomes unstable when there are equilibrium positions 
infinitesimally near to the straight equilibrium position under 
the same axial load. The correct answer to this question is 
given by replacing, in Euler’s equation, Young’s modulus by the 
so-called reduced modulus. Mr. Shanley’s analysis represents 
a generalization of the question. His procedure can be formu- 
lated as follows: What is the smallest value of the axial load at 
which a bifurcation of the equilibrium positions can occur, re- 
gardless of whether or not the transition to the bent position 
requires an increase of the axial load? The answer to this ques- 
tion is that the first equilibrium bifurcation from the straight 
equilibrium configuration occurs at a load given by the Euler 
formula when the Young’s modulus is replaced by the tangent 
modulus. In fact, one can construct sequences of equilbrium 
positions starting from any load between the two limiting values 
corresponding to the tangent and the reduced moduli. The 
inclination of the equilibrium lines representing the load as a 
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function of the deflection is steepest for the line starting from the 
lower limiting load and becomies zero for the line starting from 
the upper limiting load. The equilibrium lines have an envelope 
that starts from the lower limiting load and—at least as long as 
the stress-strain curve can be considered straight and the de- 
flection small—approaches asymptotically the load computed 
with the reduced modulus. 

“Two aspects of the question are worthy of mention: 

“(a) My original analysis, and also Engesser’s, is a generali- 
zation of the reasoning used in the theory of elastic buckling. 
Why does this not cover all possible equilibrium positions in 
the inelastic case? Obviously, it is not because of the non- 
linearity of the stress-strain relation in the inelastic range but 
because of the nonreversible character of the process. There 
are infinite values of permanent strain which may correspond to 
the same stress, corresponding to different history of the loading 
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and unloading procedure. Hence, the definition of the stability 
limit must be revised for nonreversible processes. This necessity 
was intuitively recognized by Mr. Shanley, which is, I believe, 
a great merit of his paper. 

“(b) Although the Euler formula with the tangent modulus 
does not, in general, give the maximum axial load to which the 
column can be subjected without large deflection, it is con. 
servative and therefore advisable to use this formula for practica] 
computation of column loads. As I have shown in my paper 
of 1909, also the load deflection curve that starts from the upper 
limiting load in general soon assumes a negative slope. Conse- 
quently, it is difficult to determine the actual peak of the axial 
load. It will certainly be between the two values that correspond 
to the tangent and the reduced moduli. These two values can 
be respectively designated as the lower and upper limits of the 
critical load.”’ 
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Increasing the Range of Axial-Flow Com- 
pressors by Use of Adjustable Stator Blades 


JOHN T. SINNETTE, JR.* 


National Advisory Committee for Aeronautics 


SUMMARY 


Atheory has been developed for resetting the blade angles of 
an axial-flow compressor to improve the performance at speeds 
and flows other than the design values and thus extend the useful 
operating range of the compressor. The theory is applicable to 
the resetting of both rotor and stator blades or to the resetting 
of the stator blades alone and is based on the adjustment of the 
blade angles to obtain lift coefficients below the stalling point of 


the blades. 
Tests were made on the N.A.C.A. eight-stage axial-flow com- 


pressor with three different blade settings: 


(1) With the design blade setting. 

(2) With the stator blades reset for 75 per cent of design 
speed and a load coefficient of 0.48. 

(3) With the stator blades reset for 75 per cent of design 
speed and a load coefficient of 0.65. ‘ 


The test results show that a substantial increase in the useful 
operating range of axial-flow compressors can be obtained by ad- 
justment of only the stator-blade angles. Considerable im- 
provement in efficiency at compressor speeds appreciably below 
the design speed was obtained for both stator-blade resettings. 
Also, a difference of 20 to 30 per cent in the peak-efficiency flows 
for the two stator-blade resettings was obtained, which at high 
compressor speeds is about seven times the maximum flow range 
for any given blade setting. 


INTRODUCTION 


. pe AXIAL-FLOW COMPRESSOR, in spite of the high 
efficiencies obtainable in modern designs,'~* has 
been somewhat limited in application because of its 
inherently narrow operating range. The useful flow 
range of a compressor at a given speed may be limited 
because df one or more of the following characteristics: 

(1) A rapid change in pressure ratio with small 
change in flow. 

(2) A marked drop in the efficiency with changes 
of flow. 

(3) Unstable operation, referred to as “surging,” 
which occurs when the flow is decreased be- 
yond a certain point. 

The useful speed range of a compressor may be 
limited by a marked drop in peak efficiency with change 
in speed or by a large difference between the peak- 
efficiency flow and the required flow at speeds appre- 
ciably different from the design speed. The explana- 
tion of those characteristics that limit the useful operat- 
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ing range of axial-flow compressors is the fact that the 
angles of attack, and hence the lift coefficients, for 
some of the blade rows are greatly affected by small 
changes in flow and speed; some of the blades are there- 
fore operating far from their optimum lift coefficients 
for speeds and flows appreciably different from the 
design values. 

An obvious remedy for these limitations on the useful 
range is to change the blade-angle settings with changes 
in operating conditions to maintain optimum lift coef- 
ficients. This method has been successfully employed 
in extending the high efficiency range of propeller tur- 
bines,®’ © and tests on axial-flow fans show that a pro- 
nounced extension of their flow range can also be ob- 
tained by the same method.’~’ For high-speed multi- 
stage compressors, however, the adjustment of the 
many rows of rotor and stator blades during operation 
involves mechanical problems of considerable dif- 
ficulty. The alternative method of adjusting only the 
stator blades, which is mechanically more feasible, 
has therefore been given particular consideration in the 
present investigation. 

A general method for calculating blade-angle settings 
for a wide range of operating conditions has been de- 
veloped at the Cleveland Laboratory of the N.A.C.A. 
The method is applicable to the resetting of both rotor 
and stator blades or to the resetting of stator (or rotor) 
blades alone and is based on the adjustment of the 
blade angles to obtain lift coefficients for which the 
blades will operate efficiently with special attention 
to the avoidance of blade stalling. The flow is calcu- 
lated by means of one-dimensional compressible-flow 
equations using an estimated polytropic exponent, the 
cross-sectional areas at the stations between each row of 
blades, and the air turning angle at the mean radius for 
each row of blades obtained from an empirical equation 
based on cascade tests. 

The calculation method was applied to the resetting 
of stator-blade angles of the N.A.C.A. eight-stage axial- 
flow compressor" ? for 75 per cent of design speed and 
for a series of load coefficients ranging from 0.28 to 
0.70. The load coefficient for peak efficiency at the 
design speed for design blade setting was 0.64. The 
results of tests on the N.A.C.A. eight-stage axial-flow 
compressor, with design blade settings and with stator 
blades reset for 75 per cent of design speed for load 
coefficients of 0.48 and 0.65, are presented herein. 
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SYMBOLS 


The symbols used in this report are defined here in alpha- 
betical order. Quantities that are a function of the radius are 


assumed 
stated. 


to refer to the arithmetic mean radius unless otherwise 
Where the subscripts are not associated with the main 


symbol in the definition, the significance of the subscripts is 
given in the list following the main symbols. 


A 


Aw 


Wm 


= cross-sectional area of passage, sq.ft. 
= local speed of sound, ft. per sec. 
= number of blades 
= blade chord, ft. 
= lift coefficient based on mean relative velocity 
= specific heat at constant pressure, B.t.u. per Ib. per 
= specific heat at constant volume, B.t.u. per Ib. per 
°F. 
= profile drag-lift ratio 
= standard acceleration of gravity, 32.174 ft. per sec. 
per sec. 
= mechanical equivalent of heat, 778 ft.lbs. per B.t.u. 
= constant in empirical relation for turning angle 
= polytropic exponent for expansion 
= polytropic exponent for compression 
= local Mach Number, V/a or W/a 
= rotor speed, revolutions per sec. 
= absolute static pressure, lbs. per sq.ft. 
= absolute total pressure, lbs. per sq.ft. 
= volume rate of flow, cu.ft. per sec. 
= load coefficient, cu.ft. per revolution 
= equivalent volume flow (volume flow corrected to 
N.A.C.A. standard sea-level temperature of 
59°F.), cu.ft. per sec. 
= radius to blade element, ft. 
= gas constant, ft.lbs. per Ib. per °R. 
= static temperature, °R. 
= total temperature, °R. 
= ratio of rotor velocity to axial component of air 
velocity, U/v 
= velocity of rotor (2xrn) at radius r, ft. per sec. 
= compressor Mach Number 
= absolute air velocity, ft. per sec. 
= axial component of air velocity, ft. per sec. 
air velocity relative to rotor, ft. per sec. 
= air weight flow, lbs. per sec. 
= ratio of absolute whirl component of air velocity at 
inlet to rotor row of blades to axial component of 
air velocity 
= ratio of change in whirl component to axial com- 
ponent of air velocity : 
ratio of absolute whirl component of air velocity 
(based on mean of inlet and outlet velocities of a 
row of blades) to the axial component of air 
velocity, w; + (Aw/2) 
angle of attack measured between blade chord and 
entering air velocity 
angle of attack at zero lift for isolated airfoil 
angle between compressor axis and absolute air 
velocity (referred to hereinafter as “‘air angle’’) 
angle between compressor axis and air velocity rela- 
tive to rotor (referred to hereinafter as “relative 
air angle’’) 
adiabatic exponent, ¢p/c, 
ratio of inlet-air total pressure P,; to standard sea- 
level pressure 
= loss ratio, power loss/power input (for a stage) 
= blade-profile efficiency for stage, 1 — ¢ 
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Fic. 1. 


NOTATION FOR A TYPICAL STAGE 


Station designations used in calculations. 


Fic. 2. Compressor setup with lagging removed showing loca- 
tion of pressure- and temperature-measusing stations, 


Subscripts 


polytropic efficiency for expansion 
polytropic efficiency for compression 
adiabatic temperature-rise efficiency, for constant 7, 


w > = T;) 


Ti/Tst., ratio of inlet-air total temperature 7; to 
N.A.C.A. standard sea-level temperature, Ts, = 
518.6°R. 

mass density of air, slugs per cu.ft. 

blade-element solidity, cB/2zr at radius r 

angle between blade chord and axis of compres 
sor 


ahead of compressor where velocity is negligible 
(Figs. 1 and 2) 

directly ahead of entrance guide vanes (Fig. 1) 

directly ahead of first row of rotor blades (Fig. 1) 

after last row of stator blades (Fig. 1) 

measuring station in outlet pipe from compressor 
(Fig. 2) 

flow average 

inlet to a row of stator blades for any stage (i = 4, 

inlet of following row of rotor blades (Fig. 1) 

an integer 

mean of inlet and exit values for row of blades 

at mean radius 

rotor 

stator 

standard 

at rotor-blade tip (casing surface) 

whirl component 
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INCREASING RANGE OF AXIAL-FLOW COMPRESSORS 


BLADE RESETTING THEORY 


Analysis of Problem 


The losses in any stage of an axial-flow compressor 
may be divided into losses resulting from blade-profile 
drag, blade-end effects, and passage-wall friction. The 
profile efficiency at any given radius depends on the 
profile drag-lift ratio and the ratio of the various veloc- 
ity components as given approximately by the follow- 
ing equation 


((1 + w?,,) (Aws/u Awe) (1) 


which is essentially the same as Eq. (12) of reference 1 
but is generalized to allow for cases where D/L and Aw 
are different in the rotor and stator blades. It is 
obvious from this equation that for high profile ef- 
ficiency it is important to keep (D/L)z and (D/L)s 
small. 

The value of D/L depends, among other things, 
upon the lift coefficient at which the blades are operat- 
ing and is generally a minimum at a lift coefficient 
slightly below the stalling point of the blades. When 
the angle of attack is increased beyond the stalling 
point, the drag-lift ratio increases rapidly with a 
corresponding decrease in profile efficiency. In order 
to obtain good efficiency, it is therefore important to 
avoid stalling the blades whenever possible. Because 
of blade-end effects and passage-wall friction, the 
overall efficiency of a stage will be less than the ef- 
ficiency based on the profile drag of the blade elements. 
Inasmuch as tip-clearance losses are proportional to 
the square of the lift coefficient, the lift coefficient 
that gives best overall efficiency may be somewhat 
less than that indicated by profile efficiency; the differ- 
ence will depend upon blade aspect ratio and blade-tip 
clearances. 

In a well-designed axial-flow compressor, all the 
blades will operate in the high-efficiency range of lift 
coefficients at the design operating conditions; at 
speeds and flows differing appreciably from the design 
values, however, most of the blades will be operating at 
unfavorable angles of attack, resulting in large losses. 
The effects of changes of flow on the angle of attack of a 
typical stage are shown in Fig. 3. In order to emphasize 
the main effects of changes of flow, the actual flow 
conditions have been somewhat simplified; the axial 
velocity and the radius to the element have been as- 
sumed constant throughout the stage. It is evident 
from Fig. 3 that a decrease in flow results in an in- 
crease in the angles of attack of both rotor and stator 
tows of blades, whereas an increase in flow produces the 
Opposite effect. Either change, if excessive, results in 


a large increase in the drag-lift ratio and a corresponding 
decrease in efficiency. This effect of changes of flow on 
angle of attack is present in all stages of a multistage 
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EFFECT OF FLOW ON 
ANGLE OF’ ATTACK 


ROTOR BLADE STATOR BLADE 


(a), BELOW NORMAL FLOW Vq aq >a), 


ROTOR BLADE STATOR BLADE 
b 
(b) NORMAL FLOW 
¥ 
ROTOR BLADE STATOR BLADE 
Ne 


(c) ABOVE NORMAL FLOW Vo 


Fic. 3. Effect of flow on angle of attack for a typical stage of an 
axial-flow compressor. 


compressor, but there is an additional effect—namely, 
that of density changes throughout the compressor— 
which makes multistage compressors much more critical 
to changes of flow fhan single-stage units. 

An increase in flow, for example, reduces the angle of 
attack and, in general, the density rise over a stage. 
As this decrease in density means that the axial velocity 
after the stage must increase more than ahead of the 
stage, the change in angle of attack and density rise 
in the second stage will be greater than in the first stage. 
This amplification process continues throughout the 
compressor, with the result that small changes in flow 
at the inlet to a compressor may produce extremely 
large changes in axial velocity and angles of attack in 
the last stages. 

The effect of compressor speed on performance is 
different in a multistage compressor than in a single- 
stage unit because the large density changes occurring 
have a marked effect on the axial-velocity distribution 
from stage to stage. The ratio of the mean axial veloc- 
ities v,+,/v; between any two stations, represented by 
the subscripts 7 and 7 + k, is determined by the con- 
tinuity equation 


Vi Pitk 


In the design of the compressor the ratio of the areas at 
any two stations is chosen to give the desired velocity 
ratios at the design operating conditions. Inasmuch 
as the density ratio depends upon the speed at which 
the compressor is operating, the velocity ratios and 
the angles of attack will be correct only for design con- 
ditions. From the foregoing discussion it is evident 
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STATOR BLADE 


ROTOR BLADE ; 
| 


' 
(A) STATCR BLADE ANGLE =Yo 


(B) STATOR BLADE ANGLE? > Vs 


Fic. 4. Effect of change of stator-blade angle on angles of at- 
tack for stator blades and for following rotor blades. 


that two important factors affecting the performance 
of a multistage compressor are airflow and compressor 
speed. Dimensional analysis shows that these factors 
are satisfactorily represented for a particular compressor 
by the load coefficient Q,/n and the compressor Mach 
Number U,/a. 


Control of Lift Coefficients by Blade Resetting 


For operation at conditions other than design, the 
blade angles can be reset to give arbitrary angles of 
attack for each row of blades at some point along the 
span of the blade—for example, at midspan; the angles 
of attack at other radii, however, cannot be arbitrarily 
chosen because the twist of the blades would have to be 
changed. The spanwise variation in lift coefficient 
will depend on the compressor design and on the extent 
of the deviations from the design blade settings. A 
general statement of the magnitude of the variations is 
_ therefore not possible, but the following factors should 
be considered in making an estimate: (1) At the com- 
pressor inlet, where the hub-tip diameter ratio is usu- 
ally small and blades are twisted considerably, ap- 
preciable changes in the radial variation of the velocity 
may occur in re-establishing equilibrium when the 
blade setting is changed; whereas (2) at the compressor 
outlet, the development of boundary layer along the 
hub and casing will have an important effect on radial 
distribution of lift coefficients, with a tendency to load 
the blades more toward the ends than in the middle. 
If conservative lift coefficients are chosen at midspan, 
lift coefficients at other radii should not be excessive 
and satisfactory efficiencies should be obtainable for 
moderate deviations from the design operating condi- 
tions. 


When only the stator blades are reset, it is no longer 
possible to choose arbitrarily the angle of attack at mid. 
span on both the rotor and the stator blades. The 
stator-blade angles must be set to keep the lift coef. 
ficients on both the rotor and the stator blades withip 
reasonable limits and, in particular, below the stalling 
points of the blades. The principle underlying the 
method of adjustment of stator blades alone is that ap 
increase in the blade angle for a given row of stator 
blades not only reduces the angle of attack, hence, the 
lift coefficient, on that row of blades but also on the 
following row of rotor blades (see Fig. 4). Conse. 
quently, stalling of either stator blades or rotor blades 
can be prevented by the adjustment of the proper row 
of stator blades: in the first case by the adjustment 
of the row of stator blades that are stalled and, in the 
second case, by the adjustment of the row of stator 
blades ahead of the row of stalled rotor blades. 


Selection of Lift Coefficients 

The lift coefficient at which stalling occurs depends 
upon many factors, such as blade section, Reynolds 
Number, solidity, blade-angle setting, and degree of 
turbulence. Experiments on airfoil cascades with re- 
tarded flow” ™ indicate that stalling occurs at appre- 
ciably lower lift coefficients than for isolated airfoils, 
but there is some evidence that the lift coefficients at 
which stalling occurs on compressor blades is higher 
than for straight cascades.!2 Because of the spanwise 
variation of lift coefficients and the lack of accurate 
information on the stalling of compressor blades, it 
would appear advisable to limit the maximum lift 
coefficient at midspan to about 0.8. 


Special Conditions at Compressor Inlet 


The conditions at the compressor inlet are differ- 
ent from those for the typical stage (for the usual design 
of compressors for high speeds) because the function of 
the entrance guide vanes differs from that of the typical 
stator blades. The entrance guide vanes induce an 
initial whirl to the air in the direction of rotation of the 
rotor, whereas the typical stator row reduces this whirl 
component of the air velocity. When the compressor 
is operating at load coefficients appreciably less than 
the design value, the angle of attack for the first rotor 
row becomes excessive. This condition can be some- 
what alleviated by increasing the angle of the guide 
vanes with respect to the axis so as to increase the whirl 
component of the air entering the first rotor row. Only 
a limited improvement can be realized by this method 
for two reasons: (1) Stalling of the guide vanes may 
result if they are turned through a sufficient angle to 
prevent stalling of the rotor blades; (2) turning of the 
guide vanes produces an appreciable shift in the radial 
distribution of axial velocity, which results in a small 
change in angle of attack at the tip of the rotor blades 
for a large change in guide-vane angle. An example of 
the effect of guide-vane settings on the angles of attack 
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TABLE 1 


Effect of Entrance-Guide-Vane Settings on the Angles of Attack and the Theoretical Lift Coefficients on First Row of Rotor Blades at 
Different Load Coefficients and Radii for the N.A.C.A. Axial-Flow Compressor 


Guide- Air Angle Angle of Attack on First Row of Rotor Theoretical Lift Coefficient on First 
Vane at Exit of Blades a; (Deg.) Row of Rotor Bladest 
Load Setting* Guide Vanest ean Mean 
Coefficient (Deg.) (Deg.) Hub radius Tip Hub radius Tip 
0.28 27 43 35.1 25.6 13.0 2.15 3.02 3.45 
0.28 37 53 32.1 25.0 13.0 1.91 2.91 3.44 
0.28 47 63 25.8 23.7 12.8 1.46 2.7 3.37 
0.38 27 43 28.1 21.7 10.4 1.62 2.43 2.72 
0.38 37 53 22.3 20.4 10.3 1.25 2.25 2.67 
0.38 47 63 y 6h | 17.4 9.6 0.48 1.88 2.49 
0.48 37 53 10.3 15.0 7.2 0.62 1.61 1.94 
0.48 47 63 — 18.6 8.7 5.5 — .42 1.01 1.62 
0.48 57 73 —90.5 —34.6 — 3.5 —2.89 —1.33 0.24 
0.65 27 43 5.5 9.0 2.3 0.41 1.03 1.05 
0.65 37 53 —14.0 3.2 0.7 —0.28 0.57 0.80 
0.70 17 33 i iaeg 9.9 1.6 0.68 1.11 0.93 
0.70 27 43 9 6.3 6 0.23 0.80 0.78 
0.7 37 53 —21.2 - 1.1 — 1.6 —0.50 0.28 0.48 


* The guide-vane setting is the angle between the tangent to the concave surface of the vane and the axis of the compressor. Be- 


cause there is no twist in these vanes, this angle is the same at all radii. 


The design setting was 27°. 


+ The air angle behind the guide vanes was experimentally determined for the design setting and was found to be approximately 43° 


for all radii. 


Because of the high solidity of the entrance guide vanes, it was assumed that for other settings the change of the air 


angle from the design value was the same as the change of the guide-vane setting—that is, K in Eq. (6) was assumed equal to 1.0. 
t The lift coefficients (based on assumption that blades do not stall) are used to indicate the favorablé or adverse conditions existing 
at the entrance to the rotor blades rather than actual lift coefficients that would be expected. 


and the lift coefficients at different radii on the first 
row of rotor blades is shown in Table 1, the values of 
which were calculated for -the N.A.C.A. eight-stage 
axial-flow compressor by the method given in Appendix 
I. 

For load coefficients appreciably below the design 
value, it is therefore not generally possible to avoid 
appreciable losses in the first row of rotor blades by the 
adjustment of the stator blades alone. The improve- 
ment in the performance of the succeeding stages and 
the shift of the surge limit by the adjustment of the 
stator blades should, however, make it possible to oper- 
ate the compressor satisfactorily at lower load coef- 
ficients than is possible with the design blade setting. 


Mach Number Limitations 
When the local relative Mach Number at the entrance 


to a row of blades appreciably exceeds a certain value, — 


referred to as the critical Mach Number, large com- 
pression-shock losses may result. For a given com- 
pressor Mach Number U,/a:, any~ change in blade- 
angle settings will alter the velocity diagrams and, 
hence, the local Mach Numbers through the compres- 
sor. When the compressor is operating near the design 
speed, the local Mach Numbers on many of the blades 
will be near the critical values, and changes in blade- 
angle settings (especially changes that make the veloc- 
ity diagrams asymmetrical) may cause the critical 
Mach Numbers on some of the blades to be exceeded. 
At high compressor speeds, changes in blade angles 
required to produce large changes in the flow may there- 
fore lead to appreciable compression-shock losses. On 


the other hand, the useful flow range for a single blade 
setting may be extremely narrow at these high compres- 
sor speeds. Consequently, the extension of useful flow 


range by using different blade settings may be large 
compared to the range with a single blade setting. At 
low compressor speeds, other factors than compression 
shock will usually determine how much extension of use- 
ful flow range is possible by blade adjustment. 


Method of Calculation 


For the determination of blade resettings for par- 
ticular operating conditions (given in terms of compres- 
sor Mach Number and load coefficient), flow condi- 
tions and blade resettings are calculated progressively 
from the inlet to the outlet of the compressor. The 
stations at which flow conditions are calculated are 
taken midway between the various rows of blades, 
as well as immediately ahead of the entrance guide 
vanes and after the last row of blades. The station 
designations used in the presentation are shown in 
Fig. 1. Unless otherwise stated, all flow conditions 
and blade angles refer to radii midway between the 
hub and casing. 

Because the entrance guide vanes have a function es- 
sentially different from that of the typical row of 
stator blades, the calculation procedure for the en- 
trance guide vanes is treated separately from that of 
the typical stage. When the stator blades alone are 
reset, it is necessary to reset the entrance guide vanes 
only when the load coefficient for which the blades are 
reset is appreciably different from the design value. 
The most desirable guide-vane setting may be deter- 
mined by calculating the lift coefficient on the first 
row of rotor blades for several guide-vane settings within 
the unstalled range of these vanes (see Appendix I). 
For a typical stage (consisting of a row of stator blades 
followed by a row of rotor blades), the calculations are 
carried out for different stator-blade-angle settings 
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1.00 


@ 


@ 
(2) 


VOLUME -FLOW RATIO, Q,/Q3 


fe) 


oO; 


Q,/o, A; cos Bs 


Determination of volume flow ahead of first rotor row 


Fic. 5. 


of blades. = 


until satisfactory lift coefficients are obtained on both 
the stator and the rotor blades. 

If both the stator and the rotor blades are to be reset, 
adjustment of the entrance guide vanes are unneces- 
sary unless it is desired to obtain a more favorable radial 
distribution of lift or more symmetrical velocity dia- 
gram in the first stage. The calculations for the rest of 
the compressor are then carried out on one row of 
blades at a time, using several blade angle settings on 
each row until a satisfactory lift coefficient is obtained 
on that row. 


Calculation Procedure for Entrance Guide Vanes 


For any given guide-vane setting, the flow conditions 
between the entrance guide vanes and the first row of 
rotor blades (station 3, Fig. 1) are calculated by the 
following procedure: The value of Q,/a; is determined 
from the compressor Mach Number U;/a; and the load 
coefficient Q,/n or from the flow corrected to standard 
sea-level temperature Q,/+/6 or W.V/6/6 by the rela- 
tions: 

Q:/a, = (U;,/a1) (Qi/m) (3) 
@ 


a, ~W7gRTu. V0 
W.R T; WR W, Ve (5) 
a1 Piv ygRT, Pov 


The leaving air angle 8; from the guide vanes may be 
determined either from direct tests on the guide vanes 
or estimated from design data or from cascade tests on 
similar vanes. If the value of 8; is known for one guide- 
vane setting y, its value §;’ for any other unstalled 
setting y’ can be determined with sufficient accuracy by 
the empirical relation based on airfoil-cascade tests. 


Bs — Bs’ = K (¥ — y’) (6) 


where K is a constant (depending on the solidity) be- 
tween 0.9 and 1.0. 
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The value of Q;/a; is then determined from the 
following implicit equation [derived in Appendix I ag 
Eq. (38)]: 

Qs 


a,A3 cos B3 


3 COS Bs 


where m is the polytropic exponent which is related 
to the polytropic expansion efficiency by the equation 


(m — 1)/m = [(y — 1)/y] "par. (8) 


The solution of Eq. (7) based on ,.:. = 0.90 is shown 
in Fig. 5 with Q,/Q; plotted against Q,/a,A3 cos 
The value of V3/a; is given by 


V3/a, = (Q1/a1A3 cos 83) (Q3/Q1) (9) 


and the axial component and whirl component are 
determined from the relations 


v3/a, = (Vs/a1) cos Bs = Qs/a1A3 (10) 
Vus/ai = (V3/a1) sin Bs (11) 


The density ratio across the guide vanes is obtained 
from the continuity relation 


ps/pr = Q1/Qs 


and pressure ratio, temperature ratio, and velocity-of- 
sound ratio are obtained from the polytropic relations 


Ps m ts m/(m—1) =) 
— = = — = ] 
pi - (*) (*) (13) 


For convenience of calculation, these four ratios are 


(12) 
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Fic. 6. Determination of pressure ratio, density ratio, 
temperature ratio, and velocity-of-sound ratio over entrance 
guide vanes. mp0. = 0.90. 
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plotted in Fig. 6 as a function of the variable 
0,/a1A3 cos for = 0.90. 

The whirl component of air velocity relative to the 
rotor can now be obtained from the relation 


= (Us/a1) — (Vus/a1) 


The relative air angle and the resultant relative air 
velocity ahead of the first row of rotor blades can be 
readily calculated from the two velocity components if 
the small radial component of velocity is neglected. 


(14) 


Calculation for Typical Stage 


The calculation procedure for the first row of rotor 
blades is essentially the same as that used in a typical 
stage and need not be presented separately. The 
special conditions encountered in the first row of rotor 
blades will be discussed after the general method for a 
typical stage is presented. A velocity diagram for a 
typical row of blades (stators) is shown in Fig. 7. The 
air enters the stator row of blades with a velocity V; 
at an angle 6; with respect to the compressor axis. 
The air in passing through the stator blades is deflected 
through an angle 8; — §;+,; and leaves the blades with a 
velocity V;+1 at an angle 6,+;. The air enters the fol- 
lowing row of rotor blades with a velocity relative to 
the rotor of W,+; at an angle ¢,+; with respect to the 
compressor axis and leaves the rotor row with a veloc- 
ity Wi+2 at an angle ¢:+2. The axial component of 
velocity v is shown in Fig. 7 as increasing from station 
to station, as is the case in the N.A.C.A. eight-stage 
axial-flow compressor. The value of U may also be 
changed slightly from station to station because of the 
change in mean radius. 

In order to obtain the conditions at the exit of a row 
of blades.from the conditions at the entrance to the 
row, it is necessary to know the angle through which 
the air is turned while passing through the blades. Cas- 
cade tests have indicated that, for solidities in the 
neighborhood of unity, a good approximation to the 
turning angle is given by the empirical rule 


(15) 
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Fic. 7. Velocity diagram for a typical blade row (stator) of an 
axial-flow compressor. 


where K is a constant slightly less than 1 and ay is the 
angle of attack for zero lift for an isolated airfoil. Fora 
solidity of unity, a value of K equal to 0.9 was found 
satisfactory. Although a value of 0.9 was used with 
satisfactory results in the resetting of the blades of the 
N.A.C.A. eight-stage axial-flow compressor, somewhat 
greater accuracy might be obtained by taking K as a 
function of the solidity increasing with an increase in 
solidity. Other empirical relations for turning angle, 
of course, may be used in place of Eq. (15), or, if suf- 
ficient cascade data for the blade-sections used in the 
compressor are available, the turning angle may be ob- 
tained from these data. 

The air angle §,+; leaving the stator row can be ex- 
pressed in terms of 8;, ao, and Ws (Fig. 7) by the use of 
Eq. (15) and the relation 


a, = Bi — vs (16) 


thus obtaining 
= (1 — K)Bi + + Kay 


The corresponding relation for the following row of 
rotor blades is 


= (1 — + + Keo 


(17) 


(18) 


With given stator-blade-angle setting and inlet conditions, the leaving air angle is calculated from Eq. (17) and 
the value of v ,+1/a; is determined by a trial-and-error solution of the following equation [derived in Appendix 


II as Eq. (50)]: 


a A; a; 


cos? Bi cos? Bit1 


(19) 


Inasmuch as v,+1/a will be nearly equal to v,/a:, a satisfactory approximation can usually be obtained by a few 
trials. The density ratio p;+1/p: is determined in carrying out the foregoing approximations and is equal to the 
expression in the braces raised to the indicated power in Eq. (19). 
velocity of sound ratios are determined from the polytropic relations 


The pressure ratio, temperature ratio, and 


(20) 
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The value of the exponent m’ can be determined from the estimated polytropic efficiency of compression by the 


relation 
m'/(m’ — 1) = [y/(v — 1)] 
and (unlike the exponent ym for expansion) is greater than y. The lift coefficient is calculated from the relation! 
Cy = (2Aw/c) cos Bm = (2/c) (tan B; — tan B; + 1) cos Bm (22) 
where 
tan Bn = (tan B, + tan B.+1)/2 (23) 


Eq. (23) neglects variations of axial velocity and effect of drag, but, because of uncertainties regarding the most 
desirable lift coefficient, this equation is considered adequate. 

The calculation procedure for the following row of rotor blades is essentially the same as for the stator blades, 
The inlet conditions relative to the rotor are determined from the conditions leaving the stator row and the relative 


air angle ¢;+2 leaving the rotor row determined from Eq. (18). 


equation [derived in Appendix II as Eq. (51)]: 


The value of v;+2/a; is then determined from the 


Vit1 


a \Awi 2 Leos? cos? 


and the lift coefficient is calculated from the equations analogous to Eqs. (22) and (23), with the air angles relative 


to the rotor replacing those relative to the stator. 


Special Conditions for First Row of Rotor Blades 


When only the stator blades are reset, conditions are 
somewhat different for the first row of rotor blades than 
for a typical stage, because, as has already been pointed 
out, it is not practical to control the lift coefficients on 
the first row of rotor blades to an appreciable extent 
by the adjustment of the entrance guide vanes. Con- 
sequently, when the stator blades are reset for a load 
coefficient considerably below the design value, the 
first row of rotor blades will be stalled and the method 
presented here will not be accurate for that row of 
blades. The actual turning angle and density rise 
will be less and the axial velocity leaving the blades 
higher than calculated. At a somewhat lower inlet air- 
flow, however, the axial velocity and the air angle 
leaving the following row of stator blades (station 5, 
Fig. 1), and, hence, the flow conditions in the subse- 
quent stages, should be substantially correct as calcu- 
lated. Somewhat higher accuracy, of course, could be 
obtained by basing the turning angle and density rise 
for the first row of rotor blades on cascade tests under 
similar stalled conditions. 


APPLICATION AND TEST OF THEORY 


Calculation of Stator-Blade Resettings for N.A.C.A. Eight- 

Stage Axial-Flow Compressor 

The preceding theory has been used to calculate 
stator-blade resettings for the N.A.C.A. eight-stage 
axial-flow compressor’ for 75 per cent of design speed 
for a number of different load coefficients ranging from 
0.28 to 0.70. The results of these calculations are 
summarized in Table 2. The blade resettings were 
calculated for load coefficients appreciably less than 


the design value of 0.64, because computations indi- 
cated that for certain practical gas-turbine applica- 
tions the compressor would have to operate at appre- 
ciably reduced load coefficients at speeds below the de- 
sign. The design theory, the theory of operation, and 
the complete details of construction of the N.A.C.A. 
eight-stage axial-flow compressor (Fig. 8) are given in 
reference 1. 

The changes of blade angles from the design settings 
are extremely large in some cases, amounting to ap- 
proximately 36° for the first four stator rows at the 
lowest load coefficient calculated. (See Table 2 ata 
load coefficient of 0.28.) Changes required are fairly 
large for all resettings (over 20° for some of the rows) 
because the speed for which the blades were reset was 
appreciably below the design speed. 


Fic. 8. N.A.C.A. eight-stage axial-flow compressor. 
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TABLE 2 


Summary of Calculations on Resetting of Stator Blades of N.A.C.A. Eight- 
Stage Axial-Flow Compressor for 75 Per Cent of Design Speed 


2g $8, 63 
= 
CEA OF so § 
3 Ss = Se 
Load coefficient, 0.28 
Entrance vanes 37.0 10 1.08 0.985 1.88 0.07... 
First rotor* 55%4 0 2.91 1.143 0.28 0.48 7.3 
First stator 84.9 35.9 0.28 1.013 0.95 0.49 7.2 
Second rotor 47.9 0 0.80 1.014 0.81 0.16 6.7 
Second stator 83.1 35.6 0.38 1.028 0.91 0.52 6.3 
Third rotor 45.3 0 0.80 1.014 0.82 0.18 5.8 
Third stator 80.6 35.3 0.42 1.035 0.89 0.53 5.4 
Fourth rotor 43.1 0 0.77 1.014 0.84 0.21 5.0 
Fourth stator 77.3 35.0 0.42 1.044 0.86 0.54 5.0 
Fifth rotor 42.3 0 0.72 1.015 0.83 0.24 4.3 
Fifth stator 73.3 30.7 0.48 1.051 0.84 0.53 4.0 
Sixth rotor 41.6 0 0.71 1.018 0.83 0.28 3.7 
Sixth stator 68.2 28.1 0.55 1.054 0.83 0.53 3.4 
Seventh rotor 40.1 0 0.71 1.020 0.84 0.32 3.1 
Seventh stator 63.2 25.5 0.53 1.057 0.81 0.53 2.9 
Eighth rotor 37.7 0 0.71 1.029 0.78 0.35 2.8 
Eighth stator 56.9 19.4 0.70 1.091 0.69 0.53 2.8 
Load coefficient, 0.38 
Entrance vanes 37.0 10.0 1.08 0.975 1.89 0.10... 
First rotor* 55.4 0 2.25 1.114 0.41 0.45 5.3 
First stator 79.2 30.2 0.32 1.015 0.94 0.46 5.2 
Second rotor 47.9 0 0.83 1.019 0.80 0.22 4.9 
Second stator 46.7 29.2 0.47. 1.035 0.87 0.50 4.6 
Third rotor 45.3 0 0.79 1.019 0.82 0.24 4.3 
Third stator 73.1 27.9 0.51 1.041 0.86 0.51 4.0 
Fourth rotor 43.1 0 0.74 1.019 0.85 0.27 3.7 
Fourth stator 67.9 25.6 0.55 1.055 0.82 0.52 3.4 
Fifth rotor 42.3 0 0.73 1.025 0.82 0.32 3.2 
Fifth stator 62.0 19.4 0.62 1.059 0.7 0.51 3.0 
Sixth rotor 41.6 0 0.71 1.029 0.82 0.36 2.8 
Sixth stator 55.1 15,0 0.69 1.058 0.79 0.51 2.6 
Seventh rotor 40.1 0 0.71 1.034 0.83 0.41 2.4 
Seventh stator 48.4 10.7 0.68 1.061 0.79 0.51 2.3 
Eighth rotor 37.7 0 0.72 1.051 0.77 0.45 2.2 
Eighth stator 45.2 7.7 0.70 1.078 0.72 0.52 2.2 
Load coefficient, 0.48 
Entrance vanes 47.0 20.0 1.37 0.926 2.60 0.12 es 
First rotor* 55.4 0 1. 041 0.74 0.33 3.0 
First stator 69.5 20.5 0.46 1.024 0.88 0.43 3.8 
Second rotor 47.9 0 0.80 1.029 0.80 0.31 3.6 
Second stator 65.7 18.2 0.62 1.044 0.82 0.48 3.4 
Third rotor 45.3 0 0.80 1.031 0.81 0.34 3.2 
Third stator 60.0 14.7 0.70 1.054 0.80 0.49 3.0 
Fourth rotor 43.1 0 0.79 1.036 0.82 0.38 2.8 
Fourth stator 52.1 9.8 0.77 1.065 0.76 0.50 2.6 
Fifth rotor 42.3 0 0.78 1.048 0.79 0.44 2.5 
Fifth stator 46.9 4.3 0.75 1.060 0.78 0.50 2.3 
Sixth rotor 41.6 0 0.68 1.045 0.83 0.47 2.2 
Sixth stator 41.1 -9 0.73 1.049 0.82 0.50 2.0 
Seventh rotor 40.1 0 0.61 1.042 0.86 0.51 1.9 
Seventh stator 83.3 — 4.4 0.71 1.050 0.83 0.52 1.8 
Eighth rotor 37.7 0 0.62 1.071 0.80 0.56 1.7 
Eighth stator 31.4 — 6.1 0.70 1.068 0.77 0.53 1.7 
Load coefficient, 0.65 
Entrance vanes 27.0 0 0.82 0.948 1.57 0.17... 
First rotor* 55.4 0 a 1.069 0.72 0.45 3.0 
First stator 55.6 6.7 0.60 1.025 0.86 0.40 2.9 
Second rotor 47.9 0 0.84 1.047 0.78 0.40 2.8 
Second stator 53.7 6.2 0.70 1.041 0.82 0.45 2.6 
Third rotor 45.3 0 0.72 1.039 0.84 0.42 2.5 
Third stator 48.5 3.2 0.70 1.042 0.83 0.48 2.3 
Fourth rotor 43.1 0 0.63 1.036 0.87 0.45 2.2 
Fourth stator 40.4 — 1.9 0.70 1.047 0.83 0.49 2.0 
Fifth rotor 42.3 0 0.57 1.043 0.87 0.51 1.9 
Fifth stator 33.1 — 9.6 0.70 1.041 0.86 0.50 1.8 
Sixth rotor 41.6 0 0.51 1.043 0.89 0.57 1.7 
Sixth stator 24.9 -—15.2 0.70 1.028 0.91 0.52 1.5 
Seventh rotor 40.1 0 0.47 1.045 0.91 0.64 1.4 
Seventh stator 15.5 -—22.2 0.70 1.025 0.93 0.56 1.3 
Eighth rotor 37.7 0 0.52 1.115 0.81 0.73 1.3 
Eighth stator 14.5 —23.1 0.70 1.053 0.86 0.57 1.3 
Load coefficient, 0.70 
Entrance vanes 27.0 0 0.82 0.939 1.58 0.18... 
First rotor 55.4 0 0.81 1.055 0.79 0.44 2.8 
First stator 52.7 3.7 0.50 1.017 0.90 0.39 2.7 
Second rotor 47.9 0 0.70 1.042 0.83 0.42 2.5 
Second stator 48.5 1.0 0.70 1.037 0.83 0.45 2.4 
Third rotor 45.3 0 0.65 1.040 0.86 0.45 2.3 
Third stator 42.8 — 2.5 0.70 1.038 0.85 0.48 2.1 
Fourth rotor 43.1 0 0.58 1.037 0.89 0.50 2.0 
Fourth stator 34.2 — 8.2 0.70 1.041 0.85 0.50 1.8 
Fifth rotor 42.3 0 0.52 1.046 0.88 0.56 8.7 
Fifth stator 26.4 —16.3 0.70 1.034 0.89 0.52 1.6 
Sixth rotor 41.6 0 0.47 1.046 0.91 0.63 1.5 
Sixth stator 18.0 —22.1 0.70 1.017 0.95 0.54 1.4 
Seventh rotor 40.1 0 0.45 1.052 0.82 0.72 1.3 
eventh stator 8.6 -—29.1 0.70 1.009 0.98 0.59 1.2 
Eighth rotor 37.7 0 0.50 1.115 0.80 0.83 1.2 
Eighth stator 8.6 -—28.9 0.70 1.045 0.89 0.60 1.2 


* The calculations are based on the assumption that the blades do not 
Stall. For the first row of rotor blades, where a high lift coefficient is some- 
umes obtained by the calculation method presented herein, the blades 
Would probably stall, making the actual lift coefficient and pressure ratio 
Considerably less than the values indicated for this row of blades. 


The most favorable flow conditions at the compressor 
inlet were obtained for the blade resettings for high- 
load coefficients. At the lower load coefficients the 
angles of attack and the “theoretical” lift coefficients 
were excessively high on the first row of rotor blades, 
which should result in stalling of these blades. Differ- 
ent entrance-guide-vane settings were tried, but the 
improvement was found to be small (Table 1), espe- 
cially at the tips of the rotor blades. The most desir- 
able entrance-guide-vane setting under these conditions 
is uncertain. Because of increased losses in the en- 
trance guide vanes resulting from large increases in 
guide-vane angles, the increase over design setting was 
limited to 20°. This large increase in guide-vane angle 
was used only for a load coefficient of 0.48, because the 
improvement in rotor-blade conditions at lower-load 
coefficients appeared to be so slight that increasing 
the guide-vane losses by large changes in the guide- 
vane settings was hardly warranted. 

Not only is the angle of attack on the first row of rotor 
blades excessively high for blade resettings for low 
load coefficients but the relative Mach Number at the 
entrance to these blades increases as the load coefficient 
decreases because of the small whirl component of 
velocity produced by the entrance guide vanes. Inas- 
much as the critical Mach Number is low for high 
angles of attack, large compression shock losses may 
also occur at the low load coefficients even when the 
compressor is operating at only 75 per cent of design 
speed. Because of the uncertainty regarding the 
losses through the first row of rotor blades, it is dif- 
ficult to estimate the overall performance that might be 
expected from the compressor with the stator blades 
reset for low load coefficients. 

The maximum practicable load coefficients for reset- 
ting stator blades appear to be determined principally 
by the relative Mach Number on the last row of rotor 
blades. At a load coefficient of 0.70 (Table 2) the 
relative Mach Number at the entrance to the last 
rotor is 0.83, which should be appreciably above the 
critical Mach Number for these blades even for the 
low lift coefficient of 0.50. 


Tests with Design Blade Settings and with Two Stator- 
Blade Resettings 


The effect on the performance of an axial-flow com- 
pressor of resetting the stator blades according to the 
theory presented was determined by comparing the 
performance of the N.A.C.A. eight-stage axial-flow 
compressor with the design blade setting and with the 
stator blades reset for 75 per cent of design speed for 
load coefficients of 0.48 and 0.65 (Table 2). The origi- 
nal blading for this compressor! was modified to im- 
prove the strength. The principal changes were that 
the rotor blades were tapered in thickness, larger fillets 
were used, and a semicircular cut was made at the root 
of the overhanging portion of the blade (Fig. 9). Modi- 
fication of the original blading as used in the present 
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tests are described in reference 2. The test setup and 
principal instrumentation are shown in Fig. 2. Tests 
were conducted over a range of compressor Mach Num- 
bers from 0.2 to 0.8; the design Mach Number of 0.968 
was therefore not reached in the present tests. Because 
of limitations of power and outlet-air temperature at 
high speeds, the tests were run with an inlet-air pres- 
sure of 10.16 in. of mercury absolute (corresponding 
to a pressure altitude of 27,000 ft.) and an inlet-air 
temperature of 0°F. At Mach Numbers of 0.2 and 0.3 
with design blade settings, it was impossible to throttle 
the inlet air to the desired pressure and still cover the 
normal flow range of the compressor because of im- 
perfect fitting of the inlet throttle. These data are not 
presented here because other tests have indicated that 
inlet conditions have a large effect on compressor per- 
formance at these speeds. The inlet throttle was modi- 
fied when the stator blades were reset so that the sub- 
sequent tests could be run at the selected inlet pres- 
sure. 

At the high speeds, the flow ranges for each com- 

pressor Mach Number were selected by starting with a 
low pressure ratio and progressively decreasing the 
flow until a definite surge occurred. At a compressor 
Mach Number of 0.8, however, the compressor was not 
surged except for the first blade resetting (load coefficient 
0.48) because of power limitation. At the lower com- 
pressor speeds, where no definite surging occurred, the 
flow was decreased as far as possible by closing the inlet 
throttle. 
_ A comparison of the performance characteristics for 
the different blade settings is shown in Fig. 10. Fig. 11 
presents a comparison of the peak adiabatic tempera- 
ture-rise efficiencies, the peak pressure ratios, and the 
load coefficients at the peak-efficiency points for the 
three blade settings as a function of the compressor 
Mach Number. 

The effect on the useful flow range is striking at the 
higher compressor speeds, as is shown in Fig. 10 for 
U,/a, = 0.8. At this compressor speed the useful flow 
range is small for a single blade setting. By using 
different stator blade settings it is possible to increase 


the useful flow range some seven-fold and still maintain’ 


high efficiency. At lower compressor speeds, the flow 
range for a given blade setting is much greater, so that 
the extension produced by different blade settings is 
less striking although still significant. 

For both stator-blade resettings, appreciable increases 
in peak efficiencies were obtained over those with de- 
sign angles at compressor Mach Numbers below that for 
which the blades were reset (Figs. 10 and 11), and the 
peak-efficiency curves are much flatter over the range 
of Mach Numbers from 0.4 to 0.8. This shift of the 
high-efficiency range to lower compressor Mach Num- 
bers was expected because the blades were reset for 75 
per cent of the design compressor Mach Number. At 
the Mach Number for which the blades were reset, 
however, the three different blade settings gave about 
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Fic. 11. Comparison of peak-performance characteristics of 


N.A.C.A. axial-flow compressor with design blade setting and 
with two stator-blade resettings for 75 per cent of design speed. 


Fic. 12. Close-up of N.A.C.A. axial-flow compressor showing 
effect of blade-angle resetting on clearances at tips of stator 
blades. Clearances increased at trailing edge, decreased at 
leading edge. 


the same efficiency. At this compressor speed the 
slight gain that might have been expected from some- 
what more favorable lift coefficients was counteracted 
by other unfavorable factors, such as: (1) rather large 
deviations from a symmetrical velocity diagram (Table 
2), (2) increase in blade-tip clearances (Fig. 12) and 
poor fairing of the blade bases with the casing resulting 
from changes of blade setting, (3) necessity for using 
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lift coefficients appreciably below the optimum on a 
row of stator blades to prevent stalling of the following 
row of rotor blades (or vice versa), and (4) poor 
flow conditions for the entrance guide vanes and first 
row of rotor blades for the low load-coefficient reset- 
ting. 

The stator-blade resetting for high load coefficient 
gave higher peak efficiencies and pressure ratios at all 
compressor Mach Numbers than the resetting for the 
low load coefficient. The somewhat poorer perform- 
ance with the resetting for low load coefficient is prob- 
ably due principally to factors (2) and (4) mentioned 
in the previous paragraph. The increases in tip clear- 
ances were much greater for the low-load-coefficient 
blade resetting, because the principal changes in blade 
settings were made on the inlet stages where the taper 
of the rotor islarge. The lift coefficient on the entrance 
guide vanes and the distribution of lift on the first row 
of rotor blades were also less favorable for the low- 
load-coefficient resetting. 


The peak pressure ratios for the stator-blade resetting 
for a load coefficient of 0.65 are appreciably higher than 
those for the design blade settings; whereas for the 
resetting for a load coefficient of 0.48, the peak pres- 
sure ratios are substantially the same as for the design 
setting over most of the Mach Number range. At 
compressor Mach Number of 0.8, peak pressure ratios 
were not reached with the design blade setting and with 
the resetting for a load coefficient of 0.65 because of 
power limitations. 


The tests with the two stator-blade resettings also 


’ show that a substantial change in the peak-efficiency 


flow at a given speed can be obtained by changing the 
stator-blade settings. With the stator-blade resetting 
for high load coefficient, the peak-efficiency flows were 
about 20 to 30 per cent greater over the range of Mach 
Numbers tested than with the resetting for low load 
coefficient. For t.th stator-blade resettings, the 
load coefficient obtained at the compressor Mach 
Number for which the blades had been reset were 
lower than calculated. The differences may be at- 
tributed to several assumptions made in the calcula- 
tion, especially the neglect of retardation of flow along 
the passage walls because of surface friction and blade- 
end losses. The test results indicate that the 
calculated flow should be multiplied by a flow 
coefficient between 0.85 and 0.90 to obtain the actual 
flow. 


CONCLUSIONS 


By the adjustment of only the stator-blade angles, 
substantial improvement of the peak efficiencies can be 
obtained at speeds below the design speed; by the use 
of different blade settings, the flows for peak efficiency 
can be varied so as to extend greatly the useful flow 
range of the compressor. 


APPENDIX I 
ENTRANCE-GUIDE-VANE CALCULATIONS 


An exact calculation of the flow at all radii through- 
out a multistage compressor, which involves problems 
of considerable difficulty, has not been attempted in 
the calculation procedure. For the entrance guide 
vanes, however, more detailed calculations of the 
flow conditions at different radii have been made be- 
cause the small value of the hub-tip diameter ratio at 
the compressor inlet makes the radial distribution es- 
pecially important here. 


Velocity Distribution After Entrance Guide Vanes 


An approximation of the radial distribution of air 
velocity and velocity components behind the guide 
vanes can be obtained if it is assumed that radial equi- 
librium has been established at this point (station 3, 
Fig. 1). Although this assumption is not strictly true 
because, in general some radial flow will exist at this 
station, the assumption should give a better approxi- 
mation than the assumption of constant axial velocity 
at all radii. The condition for radial equilibrium is 


Op/Or = pV,,*/r (25) 


For the present calculation, losses through the entrance 
guide vanes will be neglected in order that Bernoulli's 
equation 

Sdp/p = K (26) 
with the same constant at all radii can be used. This 


equation can then be differentiated with respect to 1, 


which gives 
V(OV/0r) + (1/p)(Op/Or) = 0 (27) 
The whirl velocity component in Eq. (25) can be 
eliminated by the relation 
V~ = V sin B (28) 


and the pressure can be eliminated by use of Eq. (27) 
to obtain 


—(1/V) (0V/Or) = sin? B/r (29) 


which, when integrated with respect to the radius, 
gives 


loge (Vi/V) = — ff" (sin®B/r)dr (30) 


The indicated integration can then be carried out if 8 
is known as a function of r. For compressors in which 
8 is constant with respect to the radius, as is approxi- 
mately the case with the N.A.C.A. eight-stage axial- 
flow compressor, Eq. (30) can be directly integrated to 
obtain 


V/V. = (r/n)~ (31) 


Eq. (30) gives the distribution of the resultant velocity 
as a function of the radius. The distribution of the 
whirl component and the axial component of velocity 
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can also be readily obtained when 8 is known as a 
function of the radius. 


Mean Velocity Behind Guide Vanes 


The continuity equation applied to the flow between 
stations 1 and 3 (Fig. 1) can be written 


pivi2ardr = p3v32ardr (32) 
where the integrations are made over the flow passages 
at these two stations. Because station 1 is chosen 
where the cross-sectional area is large and the velocity 
negligible, and the fluid is assumed to be homogeneous, 
the density at station 1 is constant, and the first inte- 
gral can be replaced by piQi. The right side of Eq. 
(32) can be expressed 


pss2ardr = JS vs2ardr = Psav. Qs (33) 


where psay. is the flow-average density. Eq. (32) then 
becomes 


Qs Qi av.) (34) 


The condition for constant total enthalpy between 
stations 1 and 3 can be written (assuming V; to be 


negligible) 


1 — (t3/t:) = V3?/2gJeph (35) 


By use of the polytropic relation between temperature 
and density, the perfect-gas laws, and the relation 
for the velocity of sound, the following density ratio is 
obtained 


{1 — [(y — 1)/2] (Vs/ar)2}”"~ (36) 


If this equation is applied at the mean radius and if the 
density at the mean radius is assumed to be equal to 
flow-average density, Eq. (36) can be substituted in 


(34), giving 


Hereinafter, the mean radius subscript mr will be omit- 
ted and all quantities will refer to the mean radius un- 
less otherwise indicated. 

If the axial velocity at the mean radius V3 cos 83 is 
assumed to be equal to’ the average axial velocity Q3/As, 
then V;/a; may be replaced by Q;/a:A;3 cos 63. If this 
substitution is made and Eq. (37) is divided by aA; 


= 


Qi @ E of (7 *) x 
a,A;3 cos B3 a,A3 cos Bs 2 


(8) 


which is the same as Eq. (7). The value of V3/a; at 
the mean radius is given by Eq. (9) and can be ob- 
tained at other radii by use of Eq. (30). 


Angles of Attack and Lift Coefficients on Firs Row of Rotor 
Blades 
The angle of attack and the lift coefficient on the 
first row of rotor blades can now be calculated as fol- 
lows: The relative whirl component of velocity at 
any radius is obtained from the absolute whirl com- 
ponent by the relation 


= (U3/a1) (Vus/a1) 


and the relative air angle ¢; is obtained from the rela- 
tion 


(39) 


tan = (Wus/a1)/(vs/a1) (40) 


The angle of attack a; (relative to the entrance veloc- 
ity W3) is given by 

as = o3 — vr 
The angle of the air leaving the first row of rotor 
blades ¢, is obtained from the relation 


= (1 — + + Keo 


and the lift coefficient is obtained from Eqs. (22) and 
(23). ‘ 


(41) 


(42) 


APPENDIX II—DERIVATION OF Egs. (19) AND (24) 


For the row of stator blades the energy equation is 


For compression with a polytropic exponent m’, 
this relation can be written 


2 2 


which, when solved for p;+1/p:, gives 


cos 83, the following relation is obtained (45) 
For the following row of rotor blades the energy equation is 
2 2 2 2 
(Wer)? _ (Wer)? _ (Ure)? _ 1) (46) 
2 2- 2 2 beta 
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where U is the rotor velocity at the radius midway between the hub and casing at the station indicated by the 
subscript. The terms containing the U’s are required because of the centrifugal-force field existing in the rotating 
frame of reference; the effect of these terms is generally small in axial-flow compressors where the change in radiys 
from station to station is small. The density ratio for the row of rotor blades may be obtained in a manner similar 
to that for the stator blades, the result being 


— 2 2 2 2 2 i/(m’—1) 
Pit1 2 aq a 
In order to use Eqs. (45) and (47) together with the relations for turning angle [Eqs. (17) and (18)] and the 
continuity relation, it is convenient to express the air velocities in terms of the axial components and the reli 


air angles. Eq. (45) for the stator blades becomes 


Pi 2 ai 


and Eq. (47) for the rotor blades becomes 


(v,/a1) (vit 1/a1) 1/(m’—1) 
cos? Bi+1 (48) 
(v4+2/a1)? 


pete _ { (2) 
Pit1 2 cos? 


aq 


cos? di+2 


The density ratio can be eliminated from these equations by use of the continuity equation, which for the 


stator blades is 


pits/pi = (Ai/A et) 


to obtain 


a a A; 2 a; cos? B; 


and 


2 /(m'— 


cos? By+1 


ay Aiti 2 Leos? 


which are Eqs. (19) and (24) of the text. 


cos? +( a ) ( ) \ (51) 
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New 


Position Reror’ ofa: Machmeter 


H. W. SIBERT* 
University of Colorado 


SUMMARY 


Asimple formula is given for the position error of a Machmeter 
in terms of the Machmeter reading corrected for instrument error, 
the altimeter position error, and the standard-atmosphere density 
and pressure ratios corresponding to the altimeter reading cor- 
rected for instrument error. 


INTRODUCTION 


URING THE WAR the U.S. Army Air Forces initiated 
the development of a Machmeter for use in 
fighter aircraft. In the future, Machmeters will prob- 
ably be used not only in fighter aircraft but also in 
heavier jet-propelled aircraft, both military and civilian. 
Since no formula for the position error of this instrument 
has yet been given, it appears desirable to derive such 
a formula. 
From the well-known relation 


= + [(y — 


connecting the total pressure p,, the static pressure p, 
the Mach number M, and the ratio of specific heats y, 
it is evident that the airplane Mach Number relative to 
the free stream can be found from the relation 


Pu = p;(1 + (1) 


(when y = 1.4), where p,, is the free-stream total pres- 
sure and ~, is the free-stream static pressure. It ap- 
pears from Eq. (1) that a Machmeter can be actuated 
by the two pressures from an ordinary air-speed head. 

The Army Air Forces Machmeter is an instrument 
with an airtight case containing two diaphragm cap- 
sules, of which one is evacuated and the other has led 
into it the pressure p,, from the opening at the end of the 
air-speed head. The pressure p, from the openings in 
the side of the air-speed head is led into the space be- 
tween the capsules and the outside case. Thus the 
movement of the diaphragms of the two capsules is 
proportional to p, and (py — pr»), respectively. By a 
suitable mechanism an A.A.F. Machmeter with a cor- 
rectly graduated dial gives values of M; according to the 
relation 


(bm — Pn)/Pn = (1 + 0.2 — 1 


when it is connected to an air-speed head. 

Stating that the dial of a Machmeter is correctly 
graduated is equivalent to saying that it has no instru- 
ment error, which is the error due to the failure of the 
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(2) 


dial readings to follow Eq. (2) exactly. Thus M; in Eq. 
(2) is either the reading of a Machmeter free of instru- 
ment error or the reading of a Machmeter corrected for 
instrument error. : 

An air-speed head when properly located with respect 
to the airplane will transmit the free-stream total pres- 
sure p,; and the free-stream static pressure p, to the 
altimeter, air-speed indicator, and Machmeter. Any 
change in the dial readings of these instruments because 
of the failure of p,, to equal p,, and p, to equal p, is 
called a position error. Instrument readings that have 
been corrected for both instrument and position error 
are called ‘‘calibrated’’ values and are designated by the 
subscript c. Thus when there is no instrument or posi- 
tion error, Pi» = Pisy Pn = Ps, M; = M, and Eq. (2) be- 
comes 


(Pu — bs)/Ps = (1 + 0.2 M,?)** — 1 (3) 


It is evident from a comparison of Eqs. (1) and (3) that 
M, as given in Eq. (3) is actually the airplane Mach 
Number. 

The determination of the position error of a given 
Machmeter attached to a given air speed-head installa- 
tion consists in finding the difference between /, and 
M, from Eqs. (2) and (3). In deriving a formula for the 
difference between M, and M; the following symbols 
will be used. 


SYMBOLS 


pressure at the opening in the end of the air-speed head 

pressure at the openings in the side of the air-speed 
head 

density 

speed of sound 

altitude 

acceleration of gravity 

p/po = density ratio 

ratio of specific heats = 1.4 in this paper 

Mach Number from Machmeter, corrected for instru- 
ment error only 

Mach Number from Machmeter, corrected for both 
instrument and position error 

= altitude from altimeter, corrected for instrument error 

only 
h. = altitude from altimeter, corrected for both instrument 
and position error 


| 


AM = M,. — M; = Machmeter position error 
Ah = h. — h; = altimeter position error 
F(M) = 3.5 log (1 + 0.2M?) 
f(h) = log p 
pi = standard-atmosphere value of p corresponding to h; ° 
pi = standard-atmosphere value of p corresponding to h; 
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Subscript 0 refers to standard-atmosphere sea-level conditions. 
Subscript s refers to free-streams conditions at altitude at which 


airplane is flying. 
Note: All logarithms are to the base e. 


DERIVATION OF FORMULAS 


High-speed wind-tunnel data have demonstrated that 
a well-designed air-speed head will transmit free-stream 
total pressure with negligible error provided the head is 
not yawed more than 15° or 20°. Thus, py, = p,, except 
for unusual airplane attitudes. 

A pressure altimeter is graduated to give pressure 
altitude h according to the standard-atmosphere rela- 
tion between altitude and pressure. Since h, is the alti- 
tude from the altimeter corrected for instrument error, 
and p,; is the standard-atmosphere value of p correspond- 
ing to h;, p; is the actual pressure transmitted from the 
air-speed head to the altimeter. In other words, p; = 
Prn.* 

From Eqs. (2) and (3) and the relations py = p,, and 


Pn = Pu 
+ 0.2M,2)35 = + 0.2M;)** (4) 


Taking logarithms of Eq. (4), 


3.5 log (1 + 0.2M,*) — 3.5 log (1 + 0.2M,?) = 
log p; — log ps 


By virtue of the definitions of F(M), f(h), AM and Ah, 
this equation can be written in the following forms: 


F(M,.) — F(M,) = f(hi) — fhe) 

* Provided, of course, that lag in the system is neglected. The 
effect of lag must be determined separately from the position 
error. However, in installations similar to those. made by the 
N.A.C.A., the effect of lag is usually negligibly small. 


F(M, + AM) — F(M,) = + Ah) — (5) 


In differential calculus it is shown that the derivatiye 
dH/dx at a point x of any function H of x is given fairly 
accurately by the formula 


(dH /dx), 2 = [H(x + Ax) — H(x)]/Ax 
provided Ax is sufficiently small. By means of this 
relation, Eq. (5) can be written as 

(dF/dM)u AM = —(df/dh),=n, Ah (6) 
provided AM and Ah are sufficiently small. From the 
definition of F(M), 

(dF/dM) = Mi = 1.4M,/(1 + 0.2M,?) (7) 

The difference in pressure dp in the atmosphere at 
two points whose elevation differs by dh is equal to the 
weight of a column of air of height dh and unit cross see. 


tion. Hence, db = —gpdh. Since f(h) = log p by defi. 
nition, it foliows from the above formula for dp that 


since a? = yp/p. Hence, for y = 1.4, 
(df/dh)n = (8) 
It follows from Eqs. (6), (7), and (8) that 
AM = [(1 + (9) 


the units of length and time in g, a, and / being the 
same. -Thus the Machmeter position error AM has been 
found in terms of the altimeter position error Ah, which 
can be obtained by the usual fly-by method. Note that 
o, and p;,/po are, respectively, the standard-atmosphere 
density and pressure ratios corresponding to h,. 
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Bending of Rotor Blade in the Plane of 


Rotation’ 


SHAO WEN YUANt 
Polytechnic Institute of Brooklyn 


SUMMARY 


A theoretical analysis of the bending of the rotor blade in the 
plane of rotation is given with consideration of the elasticity of 
the blade. The analysis also includes the bending of a blade 
with variable mass distribution along the span. The results 
lead to matrix equations that are convenient for arithmetic 
calculation. 

The application of this theoretical analysis is illustrated by 
numerical calculations. The results found in the case of a 
typical rotor blade show that the maximum deflection in the 
plane of rotation is about one-tenth of that in the flapping plane 
and that the maximum bending moment in the plane of rotation 
is about one-fourth of that in the flapping plane. When a blade 
with variable section is compared with one with an equivalent 
constant blade section, the bending moment for the variable 
section is found to be appreciably greater than that for the 
equivalent constant section in the region of the inner half span. 


INTRODUCTION 


a. BENDING OF THE ROTOR BLADE in the plane of 
rotation is of particular importance in the design 
of the present steel spar and fabric-covered blade 
since the bending rigidity in the plane of rotation is 
about equal to that in the flapping plane. The de- 
flection of the rotor blade in the plane of rotation gives, 
as observed in flight tests, a certain amount of pitching 
moment that has an adverse effect upon control. It 
is expected that this pitching moment can be reduced 
to a minimum by initially bending the blade forward 
if the amount of deflection in the plane of rotation is 
known. 

In the solution of the problem of the bending of the 
rotor blade in the plane of rotation, the differential 
equation of motion of the blade is derived from a 
consideration of the forces that act on a blade element 
of length dr. They are the so-called aerodynamic 
torque, the component of the centrifugal force per- 
pendicular to dr, the inertia force, and the Coriolis 
force. 

Since the weight of the blade is seldom distributed 
uniformly over the span, a method of calculating the 
bending of the blade with variable mass distribution 
is developed. If the mass distribution and the bending 
rigidity distribution are known for a given rotor blade 
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of variable mass distribution, the bending moments 


and the deflections can be calculated in. a manner 
similar to that used in the case of the uniformly dis- 
tributed mass. The differential equations are then 
solved by the method of collocation. The same degree 
of accuracy is obtained as in the case of constant mass 
distribution along the blade span. 

In this paper the analysis is carried out for the 
nontwisted blade with constant chord. However, for a 
twisted blade with tapered chord, only the blade 
loading expression need be modified. The expressions 
presented were developed on the basis of the first 
three terms of the Fourier series expansions because 
the consideration of more terms would lead to entirely 
unwieldy equations. On the other hand, the numerical 
example was worked out using five terms of the Fourier 
series. The maximum bending moment found when 
five terms were considered was about 20 per cent higher 
than that obtained when three terms were used. 

The results found in the case of a typical rotor blade 
show that the maximum deflection in the plane of rota- 
tion is about one-tenth of that in the flapping plane 
and that the maximum bending moment in the plane 
of rotation is about one-fourth of that in the flapping 
plane. 

When a blade with variable section is compared 
with one with an equivalent constant blade section, 
in the region of the inner half span the bending moment 
for the variable section is found to be appreciably 
greater than that for the equivalent constant section. 

The reason for this phenomenon is that the variable 
blade section has a greater bending rigidity in the 
inner portion of the blade span. It is somewhat un- 
conservative to assume an equivalent constant section, 
as done in present practice, since the maximum bending 
moment of the rotor blade in the plane of rotation 
generally occurs at about 20 per cent of the span. 


SYMBOLS 
a = slope of life curve of blade element 
Cc = chord of blade element 
EI = bending rigidity 
m = mass of blade per unit length 
Mrg = bending moment of rigid blade 
Mr, = bending moment of blade with consideration of blade 
elasticity 
r = radius to the blade element 
R = tip radius 
R, = distance between center of rotation and lag hinge 


60 = blade pitch angle 
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axial flow coefficient 

coefficient of advance 

mass density of air 

angle of rotation of blade measured from neutral 
position back 

angular velocity of rotor 


FUNDAMENTAL ASSUMPTIONS 


In the derivation of the differential equation of the 
bending of the blade with consideration of the blade 
elasticity, the following simplifying assumptions are 
made and the following definitions are used: 

(1) The chord is constant along the span of the 
blade. 

(2) The twisting of the blade is neglected. 

(3) The distance 7 from the center of rotation to 
an element of the blade does not change when the blade 
is bent. 

(4) The flapping motion is represented by the first 
three terms of a Fourier expansion: 


B —acosy — hsiny 


(5) The lag motion is also represented by the first 
three terms of a Fourier expansion: 


¢ =A, + Aicosy + Bi siny 


(6) The instantaneous lag angle of an element 
of length dr is defined as the angle ¢’ = ¢ + (2/r) 
formed by the ray from the lag hinge to the element in 
question with the reference line that connects the 
center of the rotor with the lag hinge. In the defining 
equation z is the elastic backward deflection of the 
element at 7 (see Fig. 4). 

{7) The direction of the tangent to the blade 
element at 7 is measured by the angle [ = ¢ + (dz/dr). 
This is shown in Fig. 4. 

(8) Squares and products of the small angles ¢, 
¢’, and I are neglected. 

(9) In the expressions obtained for the aero- 
dynamic forces only those corresponding to the first 
three terms of a Fourier expansion are considered. 

(10) ‘The tip loss of the blade is disregarded. 

(11) The flow through the rotor disc is considered 
uniform. 

(12) The deflection of the blade is assumed to be 
of the form Z = Z, + Z, cos y + Z; sin y. 

(13) The profile drag coefficient of a blade element 
is approximated by the following function of the 
operating angle of attack of an element: 


Cu = 52a,” 
where 6) = 0.0075, 6. = 0.8. 


BENDING OF AN INFINITELY RIGID BLADE 


When investigating the bending of a blade with con- 
sideration of blade elasticity, it is convenient to discuss 
first the bending of an infinitely rigid blade. 


Cy d) dr 
dr 
Fic. 1. Forces acting on an element of a rigid blade. 


cosy + WR 


Fic. 3. Velocities at the blade: (a) in plane of disc; (b) in plane 
containing rotor axis and blade. 


The forces acting on an infinitesimal element of 
length dr are shown in Fig. 1. The aerodynamic 
forces may be reduced to a resultant acting through 
the centerline of the span and perpendicular to the 
blade element. The centrifugal force introduces a 


term whose component perpendicular to the element 
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dr is mdrw?Ryf, and the Coriolis force is a term equal 
to 2mwBBrdr. The shearing force, the bending mo- 
ment, and the longitudinal tension along the blade 
element are represented by S, M, and P, respec- 


tively. 

The acceleration of the element of length dr normal 
to the blade and in the plane of rotation is given by 
the equation of motion: : 
mdrr¢é = dS + — Cre)dr — 

mdrw*Ryt — 2mdrrwBB (1) 


The bending moment at any point r is given 
by 


Mre = — SK Sx dsdr = — Sk Sk mrédrdr — 
Sk Sk mRyw%drdr — Su 2mwBBrdrdr + 


which represents the aerodynamic torque, the following 
assumptions are employed: 
(14) The influence of the radial velocity on the 


blade forces is negligible. 2 

(15) The acute angle ¢ between the resultant U 
of the velocities UT and UP and the rotor disc is small. 
Consequently, the following expressions can be used: 


Up ad U, Up = gU 
(16) The lift coefficient C, of the blade is propor- 
tional to the angle of attack: 
= ada, = + ¢) 


Hence, the expressions for the components of the re- 
sultant flow velocities for a blade section at distance r 


are (see Fig. 3): 


1/spC Se Sn U*(Ca, — C,)drdr (2) Ur = wr + poR sin 
In the evaluation of the last integral of Eq. (2), Up = wR — rB — pwRB cos 


Substitution of the expressions for ¢, ¢, 8, 8, U, C, and C,, in Eq. (2) and integration yield the following expres- 
sion for the bending moment of the infinitely rigid blade in the plane of rotation: 


\ 


Myo = —mu*RyRyAy + + + + (2b — + — a) + 
RyRa) + — + (Re + + » 
( — 2uR*R,ddy + uRRearbs} (82 — a)]} cos ¥ + {mw*[By(Re — RyRy) — 


+ 2 + + (25 — — (Re + (be a) 


where 
-¢ 1 \® Ifr 1 Mfr 1 


The coefficients Ao, A:, and B, are determined from the condition Mpgg = 0 atr = Ry. They are: 


2 
Ay = + Rov += + (25 — a)Oo(RRavd) + (62 — a)| + 
2mRyRay 
2 3 b 2 
+ Aer + b,*) + R*Ray (a. + 4 + ARRyvadh |} 
C 
A, = p a, | + (Rev + Yor | + a) x 


3 2R 
[ + wWR?Ray (- au) 2uR*Rayddo + | + BV -) 
2 RyRav Rey 


B, + 5626”) [2uRRav] + (252 a)6o (Rev + 
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1/Ry\? Ry | ‘| E : 


BENDING OF A BLADE WITH CONSIDERATION OF ITS 
ELASTICITY 


The forces acting on an element dr of an elastic blade 
are shown in Fig. 4. The equations of motion of the 
blade element are as follows: 

(a) Normal to the element in the plane of rotation: 


mdrr<' = '/spCU*(Cy, — Cre)dr cos — — 
2mwBBrdr cos — ¢’) — mdrrw* sin (T — + 
+ dS’ (5) 
The component P’dI perpendicular to the element 
appears on account of the curvature of the element. 
It may be noted that the equation of motion of an 
infinitely rt blade can be obtained from Eq. (5) by 


putting ¢’ = ¢. 
(d) piel to the element in the plane of rotation: 


0 = mdrrw? cos (T — fe) — 
'/spCU*(Cy — Cry)dr sin — £’) + 
2mwBBrdr sin — + dP’ (6) 


Neglecting second order terms in ¢, ¢’, and T° and 
integrating, one obtains 
P’ = '/2mw(R? — r*) (7) 
When Eq. (1) is subtracted from Eq. (5) and use is 
made of the connection ¢’ = ¢ + (%/r): 
mdrr(z/r) = — Cry)dr — 2mwBBrdr| X 
[cos (TC — ¢’) — 1] + P’dl + dS’ — dS — 
mrw*dr[sin — fc) — sin (§ — (8) 


r 
Eq -C, ¢) drcos 
drcos 


Forces acting on an element of elastic blade in plane of 
rotation. 


Fic. 4. 


With substitution of the above well-known relation 
in Eq. (9) and the new symbols X = (r — Ry)/(R - 
Ry), Z = (z/R), and k = (mR‘w*)/(2EI), the differential 
equation of bending becomes 


With second order terms in ¢, ¢’, and I omitted from 
Eq. (8) and upon substitution of the expression given 
for P’ in Eq. (7), 


dr R R dX? 
From the theory of strength of materials the following + R R *, R + 
relation between shear force and deflection is known: ok / _ RX — R,/R) 
S’ = dM'/dr = —EI(d*z/dr*) di EI (=) (10) 


The term on the right-hand side of Eq. (10) represents a rigid blade loading that is given in Eq. (1) and can 
be expressed as follows: 


Ry\‘dS 
[Ji,1 X? + Sie X + Sis] + + Jo,2X + Je,3] cos + X24 


where 
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Ry\* . Ry Ry Ry\ . Ry Ry\? 


n= 1, 2,3; and 


( R 

jue = (1 [2b — a)%r — — a)paod;] 
3 2 


2 2 
[> + pra, + > + + 4] 


- R Ry 
; = El (1 (5 
R3 Ry 2 a) 


= (R°/ED[1 — (Rv/R)}*[(0/ 


R 
+ 5 a”) |+ mRyw*(B, — 
3 
js = (1 (5 (2h a) — (52 — — | 


The coefficients Ao, Ai, and B, are given in Eq. (4). 


1 — + +  — 0/21 (018 + 


[— (26: — @)Oouao + (52 — a@)(2Ab, + + mRyw*(A; + 2aoh:)] 
— (2 — a) (3 + Pua) | 
)CR*w*)[— — 

(5 CR*w? + 529?) + (52 — a) X 


> (12) 


It must be borne in mind that the moment of this loading about the hinge should be zero for all values of y. 


This requires: 


So Se + X + Si, s)dXdX = 0, So’ (Jo,1 X* + X + Jo, s)dXdX = 0, 


So Se’ (Ja, 1 X* + + Js, 2)dXdX = 0 (13) 


The above conditions can easily be verified by substituting the corresponding values of Ji, 1, Ji,2, Ji, s.... 
This serves as a valuable check of the arithmetical calculations. 


SOLUTION OF THE DIFFERENTIAL EQUATION 


In solving Eq. (10) the form Z = Z; + Z: cos y + 
Z; sin y is assumed, where Z;, Z2, Z; are functions of 
X only. This is the same harmonic form as that of the 
loading function dS/dr. By substituting Z and the 
right-hand side member of Eq. (11) in Eq. (10), the 
following differential equations are obtained. 


a‘Z, Ry 
e(1 R(X) dX? + 
Ry\3 dZ 
Ive X J, 
az Ry\? (14) 
_ 
axe * x) ® 
R 
) 


+ + Jas (14) 


R(X) = [1 - (Ry/R)]X + (Ry/R) 


The hinged root end of the blade and the free tip 
end give the following boundary conditions: 


Z = 0, @Z/dX* = 0 (for X = 0) 
7 /dX* = 0, d*Z/dX* = 0 (for X = 1) 


| 
ial 
dX R dX? 
4) dZ; Ry 4 
2k{ 1 — —) R(X) — 
( R dX R 
where a 
_ Ry Ry Ry Ry ; 
R(X) = —(2 2—11- xX - 
R} R R R 
Ry 2 
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The second and third differential coefficients are, of 
course, proportional to the bending moment and shear- 
ing force, respectively. 


Eqs. (14) can be solved by the method of Frobenius. 
This, however, is extremely laborious because of the 
rather slow convergence of the series. A more con- 
venient method is provided by the “collocation’ 
approximations of Frazer and Jones. The assumed 


form of approximation is: 


LZ = CiX + C? x 
j=3 
Uj-) yy 
x 


(j + + 2) 


xi ( 15) 


Hence, all the independent polynominals satisfy the 


boundary conditions. J. A. Beavan found that the 
approximation was sufficiently close to the solution 
calculated by the method of Frobenius when s was 
taken as 5 so that the differential equation was satisfied 
at x = 0, 1/4, 1/2, */s, and 1. 


Substitution of the right-hand side member of Eq. (15) in Eqs. (14) yields the following matrix equations: 


Zi = + —X* + 0.3X%) + — 1.2X8 + 0.4X9 + CG! (x: 4 x") 


A;,1 A1,2 Ais Ais Ci! 
A3,1 As,2 Ass Asz,4 Ass Ci 
Ag1 As,2 Ass Ass Aas Cs 
As,1 As,2 As,s Cs! 


15 
Ji,3 
0.0625/1,1 + 0.25/1,2 + Sh, 
= 10.2 Jur: +0.5S,2 + Jus (16) 
0.5625S;, 1 + 0 75Si, 2 Ji, 3 
+ Ji, 2 + Ji, 3 


Ayi =k (1 - + 0.5 *), As,1 = 2k (1 *), Ais = —24 

Az. = —{15 + k[1 — (Ry/R)]?[0.7256 — 0.4473(Ry/R) — 0.2783(Ry/R)*}} 

As. = —{6+ k[1 — (Ry/R)]*[0.2188 — 0.375(Ry/R) + 0.1562(Ry/R)*}} 

Ags = {3 (1 Rr)" —0.0888 + 0.3691 + 0.2197 [12 +k (1 
= 24, = -{3 +k (1 - 0.3748 — 0.1996“ — 0.1753 

Ass = —{12 + k[1 — (Ry/R)]?[0.3625 — 0.375(Ry/R) + 0.0125(Ry/R)*]} 

Aus = —{3 + — (Ry/R)]?[0.3535 + 0.2005(Ry/R) + 0.153(Ry/R)*]} 


Ass = {24 + 0.8k[1 — (Ry/R)]}*} 


Ry 7 Ry Ry 3 
Ao, 4 = {6.25 = k (1 “4 0.1585 0.0784 0.0801 (*) } (Ai, 4 0) 
As,4 = —{10 + k[1 — (Ry/R)]?[0.3542 — 0.3125(Ry/R) — 0.0417(Ry/R)*}} 
Aga = —{11.25 + R[1 — (Ry/R)]*[—0.1837 + 0.0785(Ry/R) + 0.1052(Ry/R)*]} 


= {40 + 0.66662[1 — (Ry/R)]*} 


Ar 5 = {7.265 —k (1 0.05998 — 0.0284 = — 0.0316 (Ai,5 = 0) 
Ass = —{3.75 + R[1 — (Ry/R)]*[0.2868 — 0.2344(Ry/R) — 0.0525(Ry/R)?]} 
Ags = —{19.02 + k[1 — (Rv/R)}*[—0.0642 — 0.0066(Ry/R) + 0.0709(Ry/R)?)} 


Ass = {60 + 0.5714k[1 — (Ry/R)]*} 


= 


= 


By 1 


Bs, 1 
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= 
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1,1 R ( R R R R 
B,, 
4 Bs, ‘ 
. 
By, : 
Bs 
Bi, 
By, 
By, 
Bs, 
By, 
al B,, 
B;, 
. 


(15) 


fy the 
at the 
lution 
S Was 
tisfied 


(16) 


BENDING OF ROTOR BLADE 291 


= Cs(X* — X* + 0.3X*) + — 1.2X* + + 5 + x’) 


10 15 3 7 
Ce? 7 sian. ye 
6 (x X +B x)+G (x 5 X +x) 


7 
B,, 1 B,, 2 B,, 3 B,, 4 B,, 5 C3? Jo, 3 
1 Bo, 2 Be, 3 4 Bz, 5 C,? 0.0625Je, 1 0.25Je, 2 + Jo, 3 
B31 Bs, 2 Bs,3 Bs, 4 Bs, C;? = 0.25 Jo1 0.5 + Je, 3 (17) 
Bui Bar Bas Bar Bas Ce? 0.5625Jo1 + 0.75J22 + 
Bs, 2 Bs, 3 Bs, 4 Bs, 5 C? Jo,1 + Jo, 2 + Jo, 3 
Zy = Cy(X* — X* + 0.3X5) + — 1.2X5 + 0.4X%) + (= 4 x") + 
10 15 3 7 
B,, 1 “Bie By, 3 B,, 4 B,, 5 C;* Js, 3 
Bz, 1 Bo, 2 Bo, 3 Bo, 4 By, 5 C3 0.0625Js, 1 + 0. 25Js, 2 + Js, 3 
Bs, 1 Bs, Bs, 3 Bs, Bs, 5 == 0.25 J3,1 + 0.5 Js,2 + Js, 3 (18) 
Bai Ba, 2 Bas Bas Bas 0.5625Js, 1 + 0.75Js, 2 + Js, 3 
Bs, 1 Bs, 2 Bs, 3 Bs, 4 Bs, 5 C;* Js, 1 + Js, .-+ Js, 3 
B,1 = —24, Boi = {15 +k (1 [0.7498 0.4957 = — 0.2541 (=) } 
B,1 = —{6+ — (Ry/R)]?[0.3626 — 0.6625 (Ry/R)+0.3 (Ry/R)?)} 


2 2 
= {3 —k (1 [ -0.2354 — 0.3377 + 0.5731 (4) = 12+ 
Ry\* Ry 
(04 + 00 


2 2 
= 24, = -{3 +k (1 0.3805 — 0.2109 — 0.1695 } 


R 
= —{12 + — (Ry/R)]*[0.425 — 0.5(Ry/R) + 0.075(Rv/R)*)} 
= —{3+ k[1 — (Ry/R)]?[—0.1477 — 0.2111(Ry/R) + 0.3588(Ry/R)*)} 
= {24+ k[1 — (Ry/R)}3[0.4 + 0.4(R,/R)]} 
Bis = 0, Bas {6.25 (1 o.t6 - — 0.0787 
= —{10 + &[1 — (Ry/R)]?[0.3824 — 0.369(Ry/R) — 0.01339(Ry/R)*]} 
= + (1 Pr)" —0.05650 0.1758 + 0.2323 } 

= 40 +k (1 (0.s80s + 0.2858 | 

Bi. = 0, Bus = {7.265 ~ 0.0603 — 0.02907 — 0.0312 
= —{3.75 + k[1 — (Ry/R)}*[0.2999 — 0.2606(Ry/R) — 0.0393(Ry/R)*]} 
B,« = —{19.02 + k[1 — (Ry/R)]?[0.01756 — 0.1702(Ry/R) + 0.1527(Ry/R)?]} 


= {60 + k[1 — (Ry/R)]§[0.3572 + 0.2142(Ry/R)]} 


2 
Bis = 0, Bes -{ 5.004 — k (1 Fr) — 0.00986 — 0.01131 } 


B;; = {2.625 — k[1 — (R/Ry)]?(0.2161 — 0.1764(R/Ry) — 0.0397(Ry/R)*]} 
= —{48.11 + k[1 — (R,/R)]?[0.072 — 0.1721(Ry/R) + 0.1001(Ry/R)?)} 


1 
i 
>. 
x) 
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Bs,s = {83.88 + k[1 — (Ry/R)]*[0.3328 + 0.1666(Ry/R)]} 


In designing a rotor blade for bending, it is necessary to assume a bending rigidity EJ of the blade in order to 
investigate the stress in the structure. The bending stress can be calculated if the applied forces, represented on 
the right-hand side of the matrix equations, and the parameter k are known for a certain helicopter or autogiro, 
The coefficients C,’, C;', .... can then be determined from the matrix equations. The corresponding bending 
moments are: 


M = [6C31X + 12(—C3! + Cy)X? + — + + 2(6C,! — + + 
20(Cs! — + 
Mm = = [6C32X + 12(—Cy? + Cy?) X? + — + 10C,?)X* + 2(6C,? — 20C,? + + 


2(10C;? — 30C.? + 21C;2)X5 + 2(15C,? — 42C;")X® + 42C,2X7) 
= + 12(—C3 + C,)X? + 2(38C,3 — 12C,3 + 10C,)X* + 2(6C, — 20C,3 + 15C,3)X4 + 

2(10C,? — 30C.* + 21C;,*)X*® + 2(15C;? — 42C,°)X* + 42C,3X7] 
Mer = Mi + M2 cosy + Ms; sin (19) 


In case the assumed EI is not satisfactory—i.e., the structure is either too weak or too strong—a new value of 
EI may be assumed and a set of new coefficients calculated. The constants k and J,, » in the matrix Eqs. (16), 
to (18) tay be changed simply by multiplying some of the terms by the factor (EJ,.w)/(ED) ig. Hence, the 
calculation is considerably simplified. 


BLADE WITH VARIABLE Mass DISTRIBUTION 


The expressions given in the previous section are 
valid only for a uniform mass, modulas of elasticity, 
and moment of inertia distribution along the span of 
the blade. When these parameters vary along the 
span of the blade, Eq. (10) may be rewritten in the 
form: 


2n'(X) El, dX? 
Ry ” 


Ry\' @?Z Ry\'dS 
— (20) 
R/ dt Elo R/ dr Fic. 5. Lag angles and deflections of rotor blade at various 
azimuth angles VY in plane of rotation calculated for typical 


where tapered spare blade. 
moR*w? m(X) ) E(X)I(X) ' 
= » &(X xX) = 
by = = = 
E'(X)I'(X) E"(X)I"(X) j 
"X)= "(X) = 
d(EI) _ a 
ET = 
The subscript 0 represents a specified reference section * 
along the span of the blade. 02 Pm > 
The term dS/dr on the right-hand side of Eq. (20) a 2 bt 
differs from the expression given for it in Eq. 
x 


(11) only in that m in the latter is changed to 
&(X)mo. 


Fic. 6. Comparison of bending moments of rigid and elastic 
blade calculated for typical tapered spare blade. 
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In solving Eq. (20) the same method can be employed 


as that used in the solution of Eq. (10). This again 
leads to matrix equations from which the coefficients 
of X in the polynomials representing Z,, Z2, and Z; 
are determined. In the solution of the matrix equa- 
tions the mass distribution function &(X), bending 
rigidity function 7(X), and its derivatives »’(X) and 
n'(X) are specified at five stations along the blade 
span—namely, at X = 0, '/s, '/2, */s, and 1. It was 
found that this approximate solution was accurate 
enough for practical applications. 


EXAMPLES 


Numerical examples were worked out in order to 
demonstrate the significance of the bending of the 
rotor blade in the plane of rotation and to show the 
differences in the bending moments calculated for a 
constant-section blade and an actual blade with variable 
mass and moment of inertia distribution. The results 
of the calculations are shown in Figs. 5 to 8. The 
numerical work was carried out according to the theory 
just presented, but on the basis of five rather than 
three terms in the Fourier expansions. 

Fig. 5 shows the lag angles and deflections of a rotor 
blade at various values of y, the angle of rotation in the 
plane of rotation. In Fig. 6 the bending moments of 
the rigid blade are compared with those of the elastic 
blade. Fig. 7 presents a comparison of the bending 
moments along an assumed constant-section blade 
with those along a typical actual blade of variable 


° ai a2 os os o7 08 10 


Fic. 7. Comparison of bending moments of actual section and 
assumed constant section. 
4 12 
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Fic. 8. Comparison of bending moments in the flapping plane 
and the plane of rotation for assumed constant section. 


section. Finally, in Fig. 8 the bending moments in 
the plane of rotation are compared with the bending 
moments in the flapping plane. 
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Letters to 


Dear Sir: 

I have read with considerable interest the paper ‘‘On the Steady 
Flow of a Gas Through a Tube with Heat Exchange or Chemical 
Reaction”’ by Paul Chambré and Chia-Chiao Lin, which appeared 
in the October, 1946, issue of the JouRNAL. I was, however, 
surprised to find several conclusions that were emphasized 
as if they were obtained by the authors, whereas, in fact, 
they were previously presented in more detailed and expressive 
form, by Szezeniowski in his paper entitled ‘‘Flow of Gas Through 
a Tube of Constant Cross-Section with Heat Exchange... .,”’ 
which appeared in the January, 1945, issue of the Canadian 
Journal of Research. 

The point is that, although Szczeniowski’s paper is cited among 
the references by Messrs. Chambré and Lin, there is not the 
slightest mention in the text of their paper (even the footnote 
number is omitted in the text), while: 

(1) The theory developed in Section 2 follows assumptions 
and lines similar to those of Szezeniowski. 

(2) Investigations given in Section 3 are, in principle, identical 
with those of Szezeniowski. 

(3) Several conclusions given by the authors in italics and con- 
cerning the velocity of flow in a tube of constatit cross section, 
the tendency of flow to reach the sonic state, the temperature 


the Editor 


change (i.e., conditions for temperature decrease while heat is 
introduced and vice versa), etc., are simply borrowed from 
Szczeniowski’s paper, with no mention made as to the origin. 

Thus, the reader of the article is led to suppose that the authors 
have for the first time discovered a few interesting aerothermody- 
namic phenomena—which is not the case. Let me cite, for in- 
stance, the following sentence from the authors’ conclusion: 
“This paper is..... It is certainly only the beginning of a 
more extensive investigation but already it has brought out 
rather interesting results. For example, it must now be kept in 
mind that the addition of heat to a flow system does not neces- 
sarily lead to an increase in temperature. Such results are 
somewhat surprising. ... . ”” Of course, they are no more sur- 
prising in October, 1946, than they were in January, 1945, when 
they were presented in Szczeniowski’s paper. 

I should like to take this opportunity of pointing out that in- 
troducing a change of mole number in computations, however 
interesting it may be theoretically, has no practical influence if 
applied to jet propulsion. It is known that this change of mole 
number is small even in theoretically correct burning mixtures; 
in burners of jet engines, moreover, a considerable air excess must 
be applied for technological reasons. The actual amount of air 
here exceeds the theoretical by not less than three times. This 
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was proved by the successful application of Szczeniowski’s theory, 
as mentioned above, to the calculation of burners for a certain 
jet-propelled helicopter. 
WiTo_p BrzozowskI 
Chief Engineer 
Jet Helicopter Corporation 


Dear Sir: 

Controlling the Eqs. (7)-(13) in the paper by Chambré and 
Lin we have found that a term dA/A is lost on the right side of 
Eqs. (10) and (11). They should be written 


(8) + dA dn 


We discovered this when treating the case “Steady Flow at 
Constant Speed with Heat Exchange.’’ Then the Eqs. (7)—(13) 
gave (with du = 0) dp = dT = 0, which is obviously abnormal. 


BENGT JACOBSSON 
Svenska Flygmotor Aktiebolaget 


Projektavdelningen, Sweden . 


Dear Sir: 

The wide interest on the subject of gaseous combustion, as 
shown by the comments and criticisms on the paper by Chambré 
and Lin, is certainly encouraging. On the other hand, the mis- 
understanding associated with some of these comments is un- 
fortunate. For instance, Mr. Brzozowski quoted that the au- 
thors described as surprising the phenomenon of a possible de- 
crease of temperature of a flowing gas upon heat addition, while the 
authors actually meant to say that this phenomenon is by no 
means surprising. (The complete last sentence, partly quoted 
by Brzozowski from the original paper, reads: ‘Such results are 
somewhat surprising if one is not aware of the inertia forces that 
must enter the picture of rapid motion of fluids.””, Mr. Brzozow- 
ski omitted the conditional clause in his quotation.) In the text 
(col. 1, p. 539), as soon as this fact is stated, the physical inter- 
pretation is givén and its consistency with the second law of 
thermodynamics is pointed out. A similar physical interpreta- 
tion is repeated in the letter of Dr. Wood to the JouRNAL. 

I should like therefore to explain a few things regarding the 
paper under discussion. The work was completed in late 1946 
when the authors were beginning to study the problem of com- 
bustion. The chief concern was to stress the fact that the aero- 
dynamic forces must be considered in the problem of combustion 
and that the study of fuels of various types is probably not as im- 
portant. The authors then went on to the investigation of the 
effect of turbulence in the flame front, and this first part of 
the paper was left untouched. It was roughly edited and sub- 
mitted when one of the authors was leaving the Jet Propulsion 
Laboratory, GALCIT. In the meantime, we read the paper 
of Professor Szczeniowski, and inserted a reference to it. Ob- 
viously, the reference number should be in the second sentence of 
the introduction, since the references are listed in the order in 
which they are cited. We wish to apologize to Professor Szezen- 
iowski for this unfortunate omission. We wish also to thank Dr. 
Wood and Mr. Jacobsson for the correction of the other obvious 
slips, which occurred under similar circumstances. 

We regret very much that we could not have given proper 
credit to Dr. Hicks. As far as we know, the report by Dr. 
Hicks was declassified and distributed at about the same time 
the paper appeared in print. (The first declassified copy was 
received here at Brown University on October 16, 1946.) It is 


clear that Dr. Hicks and Professor Szczeniowski reached similar 
results independently, although the declassification took place a 
year and a half after Professor Szczeniowski’s paper had been 
published. There might conceivably be others who haye 
reached similar results during these years of restricted publica- 
tion. 

We wish, however, to mention the following points to explain 
why Professor Szczeniowski’s paper was not discussed in detail. 
In doing so, we hope also to make it clear that the main problems 
we have in mind are not the ones already considered. 

(1) Our main interest is in the mechanism of combustion, 
which implies other considerations besides heat transfer. Only 
the part discussed in Section 3 is connected with Professor 
Szczeniowski’s results. Even in this case, our discussions car 
be easily extended to include the case where y is not a constant, 
while his discussions are essentially based on a constant value 
of 

(2) Indeed, after the general equations [Eqs. (7)—(13)] have 
been written down, many more conclusions can be obtained in 
addition to those given in the specidl case treated in Section 3 to 
which the main criticism was directed. Some of these are given 
in Section 4, which are related to condensation and detonation 
shocks. Even in the special case of constant mole number, the 
change of cross section has interesting consequences. For ex- 
ample, in an expanding tube, the flow may depart more and more 
from the sonic state while it is being heated. In a contracting tube, 
the temperature of the gas may decrease during heating for a certain 
range of Mach Numbers less than 1/+/y. These conclusions follow 
immediately from our general equations, but do not seem to be in 
the scope of Professor Szczeniowski’s paper. 

(3) The second paragraph of Professor Szczeniowski’s Ab- 
stract puzzles us very much. It states: ‘‘The conclusion is 
reached that the velocity of flow in the tube, if the tube is heated 
or cooled, shows a natural and permanent tendency to reach the 
velocity of sound, after which the heat exchange is no longer 
possible.’’ As far as we can see, cooling does not lead the stream 
to the sonic state. 

(4) It seems that Professor Szczeniowski hesitated to con- 
clude the physical reality of decreasing temperature associated 
with heat addition. (Cf. the discussions after p. 8 of his paper.) 
So far as heat transfer is concerned, the frictional effect discussed 
by Dr. Wood is enlightening. So far as combustion is concerned, 
there is no doubt that such a phenomenon does occur. 

Although our chief aim is to study the phenomenon of com- 
bustion without immediate concern toward application, I wish 
also to say one word about a remark of Mr. Brzozowski regarding 
the practical application of our results to jet propulsion. The 
use of monopropellants in rocket tubes is certainly an example of 
jet propulsion where the change of mole number should be con- 
sidered. 

Finally, I wish to reiterate that the purpose of the paper is to 
stress the point that in combustion problems it is important to 
bear in mind the dynamical forces. The individual conclusions 
are cited to exemplify this point. One should really not be too 
concerned with them. Also, the method of approach used is 
elementary and has been used by various authors.* We would 
like to admit the priority of Professor Szczeniowski and Drs. 
Hicks and Wood concerning the two results they obtained, but 
we believe that we have considered a much more general prob- 
lem and added a few more conclusions. We still feel, however, 
that we are not very far in the problems of combustion, and we are 
humbly looking forward to some further developments in these 
directions instead of concentrating our attention on a few simple 
results already obtained. 

C. C. Lin 
Brown University 


* E.g., Bailey, N. P., JoURNAUL OF THE Agronatticat Scrences, Vol. 11, 
p. 231, 1944. 
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SUMMARY 
ustion, 

Only A simple choking device for use in transonic flow experiments 
ofessor isdescribed. It is shown that most observed shock-wave oscilla- 
nS car tions in wind tunnels are due to speed fluctuation and can be 
nstant, eliminated by using a choking device. 

value 
| INTRODUCTION 
ned in CHLIEREN OBSERVATIONS of shock waves in transonic 
“2 bs flow, such as on airfoil sections, in ducts, etc., show 
nation that the waves are generally unsteady, oscillating over 
er, the appreciable distances. This apparent unsteadiness of 
or ex- shock waves has often been accepted as a typical 
i more phenomenon of transonic flow. This view was sup- 
a ported by some flight observations of severe buffeting 
follow in transonic flow. 

» be in During investigations of the supersonic zone on air- 
foil sections at the GALCIT, it was realized that, in 
s Ab- order to increase the accuracy of measurements, an 
was: attempt should be made to stabilize the flow. It was 
ch the then found that most observed shock-wave oscillations are 
onger introduced by the wind tunnel and can be stopped by 
‘ream means of a “‘choking’’ device downstream from the 
model in the test section or diffuser. A simple device 
oa of this sort has been used with excellent results for 
per.) some time in the 2- by 20-in. Transonic Tunnel at the 
ussed GALCIT. The present paper describes this device 
rned, and its application. 
The stabilization of shock waves in wind tunnels 
a and ducts is of considerable practical and theoretical 
‘ding importance. The usual rapid motion of a shock wave 

The causes inherent errors in pressure distribution and force 

le of measurements because of the difficulty of obtaining 

a these measurements with an instrument of small time 

is i lag. This shock-wave instability also introduces 

te nonstationary effects that are difficult to evaluate. 
sions There are additional important consequences in con- 

* too nection with the comparison with free-flight data and 

d is especially shock-buffeting problems. From the ana- 

a lytical point of view, it is, of course, important to know 

but whether the positions of the shock waves in flow with 

rob- local supersonic regions are well defined or not. 

ver, The large-scale oscillations of the shock wave that 

a were observed in the 2- by 20-in. Transonic Tunnel - 

ai were, in general, due to fluctuations of the tunnel flow 

(.e., to fluctuations in mean speed, not to turbulence). 

P The shock oscillations appeared to be less pronounced 

sity 


Received November 6, 1946. 
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California Institute of Technology 


ions in Wind Tunnels 


than had been noted in previous observations else- 
where. This was attributed to the comparatively 
steady flow in the tunnel and led to an attempt to 
decrease the fluctuations and pulsations even more. 

To obtain constant mass flow in the tunnel, a con- 
striction of the test-section area was established down- 
stream of the model, and the tunnel operated in such a 
fashion that this constriction ‘‘choked,” i.e., sonic 
velocity was attained there. The Mach Number in 
the test section is then, as in supersonic wind tunnels, 
determined by the area ratio of test section to throat. 
This procedure proved successful and the shock waves 
became steady. This second throat also provides an 
excellent method of speed control. A change in the 
throat area can be calibrated directly in Mach Number 
regardless of Reynolds Number; or, at a given throat 
area, the Reynolds Number can be altered by changing 
the pressure in the tunnel while the Mach Number 
remains fixed by the area ratio. 

The authors believe that the establishment of such 
a throat area is necessary for all measurements in the 
transonic range, except in cases in which the unsteadi- 
ness of flow and the subsequent loss of accuracy in the 
measurements are of minor importance. The im- 
portant parameter in the transonic range is not the 


Fic. 1. Sketch of test section showing the choke, composed ot 
flexible airfoil (1), with rotating center piece (2), downstream 
from the model (3). The position of the center piece is shown 
on the scale (4). 


295 


° 
Place a 
id been 
D have 
ublica- 
ES 
ag 
= 
2 
4 
| 


296 JOURNAL OF THE AERONAUTICAL SCIENCES—MAY, 1947 


Fic. 2. Test section of the GALCIT 2- by 20-in. transonic 
wind tunnel, showing the motor drive for the choking unit speed 
control. i 


Mach Number, /, but the parameter 8 = V1 — M?. 
The use of 8 instead of M shows immediately that the 
transonic speed range is not so small as it appears 
when the Mach Number only is considered. The 
necessary precision of Mach Number setting increases 
as M approaches unity, since 

dM 


APPARATUS 


Choke 


The variable choking device is shown in Fig. 1. 
Fig. 2 is a photograph of the test section of the 2- by 
20-in. Transonic Tunnel at GALCIT showing the 
motor drive for the choking unit. Using this device, 
the Mach Number in the test section can be varied 
continuously from 0.700 to 0.905, and values of the 
Mach Number may be reproduced with an accuracy 
of 0.5 per cent. 

The choking device is a flexible airfoil formed of two 
pieces of thin spring steel, hinged at both ends, and a 
rotating center piece. The leading-edge hinge pin 
passes through the tunnel walls and is fastened securely 
to them, while the trailing-edge pin is allowed to slip 
between the tunnel walls. As the center piece is 
rotated, the maximum thickness of the choke is varied, 
changing the throat area of the tunnel and, hence, the 
speed in the test section. 

It should be noted that the choke need not be sym- 
metrical. In fact, the first experiments leading to the 
present design were made by spanning the tunnel with 
a streamlined ‘strut section and merely rotating the 
strut. Observation with the tunnel schlieren apparatus 
does not reveal any asymmetrical flow conditions in 
the test section with such an asymmetrical choke. 


Mach Meter 


A simple and inexpensive form of Mach Number 
meter (Fig. 3) has been developed in conjunction with 
the speed control. The meter consists of a large board 
on which are mounted two mercury tubes, one of 
which is a barometer-type tube, measuring absolute 
pressure, while the other is a U-tube formed by a 
capillary tube and a large diameter reservoir. The 
static pressure is measured by the barometer tube, and 
the difference (bo — p) between stagnation pressure 
(po) and static pressure (p) is measured by the U-tube. 
The scales of these two tubes are also the scales of a 
nomograph, the slanting scale of which is also mounted 
on the manometer board. Movable pointers run along 
the scales and the indicating line of the nomograph (in 
the form of a fine thread) is attached to the pointers, 
The Mach Number is read directly from the slanting 
scale, which has been calculated* through use of the 
well-known isentropic pressure relation: 


po/p = {1 + - 


* Calculation of the nomograph and the design of the meter 
are due to P. Johnson. 
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Fic. 3. The nomographic Mach Number meter. 
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Fic. 4b. Schlieren ph —- of flow past 3-in. by 12 per cent circular arc airfoil at Ma 
Knife edge normal to mean flow direction. 


1/100 secs. 


The Mach Number scale is such that variations in the 
third decimal of the Mach Number are easily dis- 
cernible. 


RESULTS 


The type of speed control described, in addition to 
providing maximum flexibility of operation, serves to 
decrease the oscillating nature of the shock wave. In 
the past, rapid fluctuations of the free-stream Mach 
Number due to diffuser separation and variation in 


Schlieren photograph of flow past 3-in. by 12 per cent circular arc airfoil at Mo 
Knife edge normal to mean flow direction. 


= 0.844. Tunnel choked. Exp., 1/100 


= 0.844. Tunnel not choked. Exp., 


other losses throughout the tunnel have caused the 
shock wave to oscillate over as much as 20 per cent of 
the airfoil chord at a given tunnel-speed setting. By 
creating a supersonic region (through use of the choke) 
at the exit of the test section, the disturbances from 
the diffuser are not allowed to travel upstream. Then, 
since the mass flow through the test section must 
remain constant by virtue of the sonic speed at the 
throat, it follows that the free-stream Mach Number 
in the test section must also remain constant and, thus, 
oscillations of the shock wave are kept to a minimum. 
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ieren photograph of flow past 3 


Schlieren photograph of flow past 3-in. by 12 per cent circular arc airfoil at Mo = 0.902. Tunnel not choked. 


Schl 


Fic. 5a. 
Fic 5b. 
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The decrease in shock-wave fluctuation is graphically 
“4 illustrated by Figs. 4 and 5, which were taken with 
+s continuous light at an exposure of '/1 sec. Figs. 4a 
s and 5a show the shock wave on the airfoil with the 
tunnel running at the choked condition, while Figs. 
od 4b and 5b show the same shock wave in the unchoked 
ee q tunnel. 
The effect of the choking speed control on pressure 
~ distribution surveys is as striking as that on shock- 
wave oscillations. Pressure-distribution measurements 
show none of the rapidly fluctuating qualities commonly 
was =< associated with transonic pressure surveys. Pressure 
surveys may be compared directly with schlieren photo- 
graphs. A typical comparison of a pressure measure- 
ment and a schlieren photograph is shown in Fig. 6. 
The large shock ‘thickness’ shown in Fig. 4a in- 
dicates that the shock wave is not completely motion- 
less. This is evident from Figs. 7 and 8, where a 
comparison is made between schlieren photographs 
taken with a spark (Figs. 7a and 8a) at an exposure of 
4 X 10 sec. and those taken with continuous light 
(Figs. 7b and 8b) at an exposure of '/9 sec. In order 
. to determine the possible cause of the remaining 
oscillations of the shock wave, a curve of shock wave 
position vs. the transonic parameter 8 = 1/1 — M? 
(Fig. 9) has been plotted. The data from which this 
curve is drawn were taken from schlieren photographs 
made with continuous light; the position of the shock 
Fic. 6. Comparison between pressure distribution measure- wave is easily read from such photographs. 
ments (top) and schlieren photograph. Knife edge normal to This curve may be used to calculate approximately 
mean flow direction. Exp., 1/100 secs. 
the change in throat area required to cause the shock 
“thickness” of Figs. 7b and 8b if the assumption is 
made that the thickness of the shock is due simply to 
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1/100 Fic.7a. Schlieren photograph of flow past 3-in. by 12 per cent circular arc airfoil at M. = 0.844. Tunnel choked. Exp.,4 X 107 
secs. Knife edge normal to mean flow direction. 
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Fic. 8b. Schlieren photograph of flow past 3-in. by 12 per cent circular arc airfoil at Mo = 0.902. Tunnel choked. Exp., 1/100 
secs. Knife edge normal to mean flow direction. 


the fact that the shock is moving during the exposure 
of the film. Using simple one-dimensional gas dy- 
namics, the variation of Mach Number with throat 
area may be determined; thus: 


(y+) 


A* \y+1 


or 
y+1 
m2 —1)/2(y—1) 
(1 +7 wr) 


dx/dM may be measured from Fig. 9; dA* may then 
be computed for any apparent shock thickness. 
Computation shows that a change of 0.001 in. in 
all dimensions at the throat could account for a Mach 
Number variation sufficient to cause the apparent 
“thickness” of the shock wave. Changes of this 
magnitude can be due to fluctuations in the thickness 
of the boundary layer. Thus, with the present ar- 
fangement it is difficult to determine how much of 
the remaining shock oscillation is ‘‘real’’ and how much 
is due to the wind tunnel. There is little doubt that 
a certain amount of the shock oscillation must be due 
to boundary-layer fluctuations on the model itself 
and, thus, is “‘real’’—that is, is independent of the 


Swoce Wave Position V3 | 
Arc Arron 


(OATA FROM PHOTOGRAPHS 
wiTe CONTINOUS FOS/TION 
or wave omy) 


—~ FER Centr of 
g 


7s 


Swoca Wave 


oo 
Free Srream —~ Me 
Fic. 9. 


tunnel characteristics. Subsequent investigations must 
show the magnitude of these oscillations. At this 
time, real shock-wave fluctuations have been definitely 
observed only when the base of the shock wave almost 
coincides with the transition point of the boundary 
layer on the model. 
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CONCLUSIONS 

The speed control of wind tunnels operating in the 
traysonic range must be considerably more precise 
than that of subsonic wind tunnels. The use of a 
choking device to stabilize and control the Mach 
Number appears to be the best choice for a transonic 
wind tunnel. It should be noted, however, that the 
use of such a device in a tunnel requires the expenditure 
of considerably more power than is ordinarily necessary. 
The energy ratio of a tunnel operated with a choking 
device is thus far lower than that of an ordinary sub- 
sonic tunnel. Such an increase in power cansumption, 
however, is well warranted by the advantages to be 
gained. If a transonic wind tunnel is to be equipped 
with such a device, it should be remembered that the 
speed fluctuations due to oscillation of the wall bound- 
ary layer can be reduced by using large test-section 
cross sections. 


The conditions prevailing in a choked transonic wind 
tunnel are, in general, nearly the same as those present 
in a normal supersonic tunnel. The fact that the 
quality of the flow downstream from the throat of a 
transonic tunnel is of no importance merely represents 
a greater latitude in the design of such a tunnel and 
makes the use of a variable throat area much simpler 
for operation in the transonic region than in the super- 
sonic region. ‘ 

The results presented here show that observed 
buffeting phenomena in transonic flight are not con- 
nected with the hitherto observed shock oscillations in 
wind tunnels. The buffeting is likely to be caused 
by a combination of structural oscillations and air- 
speed fluctuations rather than by any inherent charac- 
teristics of transonic flow. Further research is evi- 
dently needed to clarify this point. 
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The Rate of Descent of Parachutes from 
Various Altitudes 


KENNETH E. PENROD,* GEORGE L. MAISON,* ann JAMES E. McDONALDt 
Boston University School of Medicine and Iowa State College 


SUMMARY 


One-hundred and eleven observations were made on the descent 
times of 200-lb. torso dummies attached to 28-ft..silk, 28-ft. nylon, 
and 24-ft. nylon parachutes. These observations were made 
from various launching altitudes up to 40,000 ft. A CaS value 
was computed for each descent from the formula CaS = 2W/pV?, 
where W = weight of the dummy plus parachute and pV? was 
calculated from the observed atmospheric conditions. The mean 
CS values thus computed were: for 28-ft. silk parachutes, 531; 
for 28-ft. nylon, 483; and for 24-ft. nylon, 336. 

By using the mean CzS value computed for each type of para- 
chute, a formula has been developed to predict the time of de- 
scent from any altitude to any altitude if the type parachute, 
weight of the man, and certain atmospheric conditions are known. 

An improved method of comparing performance of experi- 
mental parachutes is presented which negates the necessity for 
attempting to duplicate a given density altitude on successive 


flights. 

A method is presented for calculating velocity at any instant 
in parachute descent. Landing velocities at different ground 
levels and with different types of parachutes and different weight 


men are discussed. 


INTRODUCTION 


KNOWLEDGE OF THE TIME REQUIRED to reach the 

ground when a parachute is opened at a specified 
altitude is sometimes highly useful. Moreover, in- 
formation concerning the rate of descent is necessary 
for the calculation of landing impact when the ground is 
reached. Heretofore, all tables of rates of parachute 
descent have been based on calculations extrapolated 
from a few observations of parachute descents from near 
the ground. It has been presumed, but never proved, 
that the rate of descent at any given altitude would vary 
inversely with the square root of the density ratio of 
the air at that altitude. In connection with tests con- 
ducted for the measurement of impact forces of para- 
chute openings at altitude, the occasion was presented 
to observe experimentally the effect of altitude on para- 
chute descent. 


METHOD 


The technique of release of the 200-lb. hard rubber 
torso dummies was described in a previous report.! 


Received January 10, 1947. 

*Departments of Physiology and Pharmacology, Boston 
University School of Medicine. Formerly at the Aero Medical 
Laboratory, Air Matériel Command, Wright Field. 

t Department of Physics, Iowa State College. 


Briefly, the parachutes were opened by 25-ft. static 
lines attached to the airplane. Parachutes were re- 
leased at approximately 7,000, 15,000, 26,000, 33,000, 
and 40,000 ft. density altitude, and at each release ac- 
curate data as to indicated altitude, air speed, and 
temperature were kept. Flight instruments had been 
previously calibrated by the Flight Test Branch of the 
Air Matériel Command, Wright Field. In addition to 
the usual aircraft altimeter, the altitude of release was 
recorded by a radio altimeter and a recording barograph. 

The timing was done visually by ground observers 
for the lower altitudes, and for the higher altitudes radio 
contact between the plane and the ground was main- 
tained and stop watches were started at a signal from 
the plane. The parachutes were released in succession 
with sufficient intervals between to allow individual 
timing. The watch was stopped with the first contact 
of the dummy with the ground. The ground was a dried 
desert lake bed, the altitude of which was +2,800 ft. 
The location was near Kingman, Ariz.; the period of 
observation was November 30, 1944, to February 10, 
1945. 


RESULTS 


Since the primary purpose of these parachute drops 
was to measure and compare opening shock forces of 
different sizes and kinds of parachutes, a given drop 
usually included a mixture of parachutes. Thus, the 
number of observations on a given type of parachute 
for any one flight varied from two to six. 

Some variation existed among altitudes of release 
with different flights. For instance, the calculated 
density altitude at release of the silk parachutes at the 
so-called “‘26,000-ft.”’ level ranged between 25,850 and 
27,050 ft. density altitude. Since there is good reason 
to presume, as will be shown later, that the rate of de- 
scent of a parachute is related to the density of the air, 
it is misleading to group all the observed descent times 
from one approximate altitude. Also atmospheric 
conditions between the altitude of release and the 
ground will vary from day to day. The observations 
herein reported were made at the same geographic 
location but over a period of some 10 weeks. There- 
fore, it should prove more fruitful to treat each flight 
individually, taking into consideration the actual ob- 
served atmospheric conditions for that day. To this 
end the Kollsman number, or altimeter setting, at take- 
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off and landing and the atmospheric temperature at 
each 2,000-ft. stratum were recorded for each flight. 

Terminal velocity of a freely falling body is reached 
when its aerodynamic drag is equal to its weight. This 
relationship is expressed thus: 


D = W = C4S(pV;7/2) (1) 
where 
D = drag (Ibs.) 
W = weight (Ibs.) 


Ca = drag coefficient 

S = drag area (sq-ft.) 

p = air density (slugs per ft.*) 

V, = terminal velocity (ft. per sec.) 


If the density of the atmosphere did not vary with 
altitude, then any falling body, as soon as it had been 
accelerated to the terminal velocity specified by Eq. (1), 
would fall with that constant velocity throughout its 
entire descent, and the problem of predicting total time 
of descent would be comparatively easy. But in the 
actual atmosphere, air density is not constant in space, 
or in time, so the terminal velocity of a body falling 
through the earth’s atmosphere becomes a variable 
depending not only on its shape and weight but also on 
the value of the air density at each successive point in 
its descent. Since a body falling through the atmos- 
phere is continuously moving into regions of greater 
density, Eq. (1) tells us that this is equivalent to saying 
that it is always falling from regions of higher terminal 
velocity into regions of lower terminal velocity; so, 
strictly speaking, we must conclude that in all but the 
first few seconds, during which the chute is opening, this 
steady deceleration implies that the actual velocity of 
descent must always be slightly greater than the theo- 
retical value determined by Eq. (1). That is, the drag 
force, D, must always be slightly greater than the 
weight, W, with the net effect that the body undergoes 
a steady deceleration, as is actually observed to occur. 
The rate of increase of density with decreasing altitude 
is, however, gradual enough to make this slight dis- 
equilibrium and its consequent velocity lag negligible 
in our calculations. Thus, in summary, Eq. (1) pre- 
. sents the following picture of the variation in rate of 
descent of any large body (if we agree to equate D and 
W, and thus to neglect any velocity lags): For a body 
of given weight, W, drag area, S, and drag coefficient, 
C4, the velocity at any point in its descent is inversely 
proportional to the square root of the density. From 
this basic relationship, formulas expressing the depend- 
ence of total time of descent on atmospheric conditions 
and on parachute type will now be developed. 

The weight of our dummies, with parachutes at- 
tached, was constant at 225 lbs., plus or minus a few 
pounds. The pV? values would be expected to vary 
with the changes in atmospheric conditions from flight 
to flight. These pV? values were computed from the 
data of each flight as follows: 


Let 
pV? = k = 2W/C.S 
Then 
V = Vk/p 
also | 


1/V = Vp/k 


If, now, we imagine the atmosphere to be divided 
into a series of strata, each 2,000 ft. in vertical thick- 
ness, then the time required to pass through the ith 
stratum will be 


t; = 2,000/V; 


The sum of these 2,000-ft. stratum times represent 
the total time 


= 2,000 2,000 i 
V; V; 


1/V = Vp/k 


But 


2,000 - 
Lt; = 2,000 Pi = 2 


Now from the General Gas Law we have 
P = pRT 
in which - 


P = pressure in inches of mercury 
p = the air density expressed in slugs per ft.* 


R = the gas constant for dry air (in suitable units) 
T = air temperature in degrees Kelvin (°K.) 
Then p = P/RT and 
Vp = VP/RT = (1/VR)VP/T = 
3.38639 X 104 X 1.951875 X (P/T) ra) (3) 
2.8886 X 10° 


0.15126 P/T 


This expression can then be substituted in Eq. (2), 
yielding 


sy, = 21000 X 2 
i T, 
(4) 


where P;/T;, represents average pressure and tempera- 
ture conditions for the ith 2,000-ft. stratum. 

If the observed time of descent for each parachute is 
now substituted in Eq. (4), we may rearrange and solve 
for k: 
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TABLE 1A 


Computed pV? and CzS Values for Each Observed 24-Ft. Nylon 
Parachute Descent 


TABLE 1B 
Computed pV? and CaS Values for Each Observed 28-Ft. Nylon 
Parachute Descent 


Observed Observed 
Flight Density Time, Flight Density Time, 

No. Altitude Sec. pV? CaS No. Altitude Sec. pV? CaS 

12 15,550 477 1.27 353.5 9 7,750 204 1.07 419.0 

i 464 1.385 334.6 204 1.07 419.0 
454 1.41 320.3 214 0.98 461.1 

474 1.29 349.1 213 0.99 456.4 

477 Be - 353.5 14 7,800 212 1.00 450.9 

456 1.39 323.0 203 1.09 413.6 

16 15,800 455 1.43 314.0 224 0.89 503.3 

472 1.33 337.8 219 0.94 481.3 

426 1.64 275.2 4 16,000 605 0.83 540.2 

468 1.35 332.4 540 1.05 430.2 

23 15,650. 465 1.33 339.1 538 1.05 426.9 

440 1.48 303.6 6 15,875 557 0.93 484.4 

504 1.13 398.2 559 0.92 488.1 

19 26,300 828 1.24 363.8 543 0.98 460.6 

785 1.38 326.3 S 15,750 552 0.99 454.5 

36 26,400 783 1.42 316.2 578 0.90 498.3 

789 1.40 321.2 563 0.95 472.7 

787 1.41 319.6 13 15,850 565 0.93 482.3 

815 1.31 342.7 612 0.80 565.3 

34 32,740 924 1.47 305.9 21 15,950 568 0.93 485.4 

925 1.47 306.8 574 0.91 495.6 

1,047 1.15 393.0 571 0.92 490.2 

40 39,000 1,150 1.22 370.1 5 27,050 ' 927 1.02 440.3 
1,149 1.22 369.5 998 0.88 510.2 

Mean 1.348 336.23 7 26,600 949 0.94 476.2 
963 0.92 490.2 

Standard Deviation 0.116 29.1 i 25,850 990 0.86 593.9 
995 0.86 520.2 

Referring again to Eq. (1), it is seen that C,S values 958 0.93 481.8 
may be then calculated by the equation 15 26,360 ae oe 
CaS = 2W/pV? (6) 968 0.90 498.3 

19 26,300 976 0.89 505.6 

where the value of pV? has been obtained from Eq. (5) 954 0.93 482.8 
and where W is the known weight of dummy plus para- 36 26,400 907 1.06 424.5 
C.S values have been calculated for 24 observed de- . — 1,460 on 565.3 
scents of 24-ft. nylon parachutes, for 46 28-ft. nylon, 1,428 0.84 540.9 
and 41 28-ft. silk parachutes. These data are pre- 18 40,000 1,410 0.86 525.7 
sented in Tables 1A, 1B, and 1C. — Sa — 
1,342 0.94 476.2 

1,386 0.89 507.9 

COMPUTATION OF APPROXIMATE DESCENT TIME 40 39,000 1,350 0.88 510.2 
Using the above relationship, it is possible to predict 
with reasonable accuracy the time it will take a man to 0.908 482.83 

from any altitude if we know his an , 

teach the ground Standard Deviation 0.077  +38.75 


weight, the type of parachute used, and the ground- 
level elevation. 
Transposing Eq. (5), 


302.522 
descent time = = (7) 


If the mean C,S values for each of the three standard 
types of parachutes are used, the formulas become: 
For 24-ft. nylon parachutes: 


302.522-/ P;/T; 


(8) 


descent time = 


V 2w/336 


For 28-ft. nylon parachutes: 
302.522-/ P./T; 


descent time = — (9) 
V 2w/483 
For 28-ft. silk parachutes: 
descent time = = (10) 
V 2w/531 


The value for 2+/ P:/T; is a function of the space 
through which the parachute travels. If the ground 


de 


4 
| 
P, 
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i a Pe to 1,000 ft., we would read the value for 20,000 ft. from = 
ce Computed pV? and C,S Values for Each Observed 28-Ft. sik the graph, 2.77 and subtract the 0.15 for the 1,000-ft. 
landing elevation, getting 2.62 as the appropriate value 
The value for shown in Fig. 1 was computed 
No. Altitude Sec. pv? CaS from U.S. Standard Atmosphere tables. Hence, the 
1 7,500 204 0.87 517.8 error of prediction should be in proportion to the Py 
‘ 210 0.99 455.0 variation of the actual atmospheric conditions of the " 
ts 233 0.80 560.4 day from the standard atmosphere. However, this 
" 9 7,750 246 0.74 608.9 error is probably less than the errors introduced by the 1000 
: on pei pone variation among parachutes. It is seen from the data ’ 
14 7,800 242 0.77 587.5 
' 213 0.99 455.5 of Table 1 that identical parachutes dropped in succes- 
f 240 0.78 577.7 sion from a plane at a constant altitude will vary con- * 
“ 4° 16,000 583 0.90 501.7 siderably in the time they take to reach the ground, 
: pc rn : = ron This variation can be partially accounted for by rising hi 
; 572 0.88 514.4 and descending air currents and the amount of pen- \ 
613 0.77 586.7 dulum swing; hence, spilling of air. For this reason it 
8 15,750 562 0.95 471.2 is possible to predict only approximately the actual 400% 
’ “ 600 0.84 537.0 time of descent of a lone parachute. Therefore, the . 
: 602 0.83 540.9 figures calculated by the above method should be taken Seo. 
i, to represent the mean, or average, descent times of 
: 558 0.96 470.2 groups of identical parachutes and loads. ee 
607 0.81 556.2 adopte 
4 552 0.98 460.1 COMPARISON OF PARACHUTE DESCENT TIMES 
21 15,950 609 0.81 557.6 
586 0.87 516.6 For the purposes of comparing various types of 
5 27,050 998 0.88 510.2 parachutes it is believed that a computation of C,§ sacly 
1,115 0.71 636.5 values would yield better results than does a simple vue 
comparison of observed descent times. Such computa- 
? - 27,600 970 0.91 497.2 tions would relieve the necessity for attempting to du- 
; 990 0.87 517.8 plicate pressure altitudes on different flights, as well as ay 
15 26,360 1,022 0.81 555.6 
1,032 0.80 564.6 a altitu 
19 26,300 984 0.88 513.7 / flights 
964 0.91 493.4 38 a ye 7 lation 
¥, 10 39,400 1,585 0.67 666.7 36F in Calculating Parachute 7 (a) 
1,372 0.90 499.5 Descent Times 
11 40,050 1,389 0.90 499.5 32 with ¢ 
: 1,388 0.90 499.5 30 / altituc 
1,389 0.90 499.5 28 the fis 
hr 1,448 0.83 542.8 = / if sucl 
18 40,000 1,417 0.85 531.3 26 7 
1,461 0.80 564.6 24 
Mean 0.855 530.74 22 
Standard Deviation +().077 +48. 63 20 7 f 
i 16 VA set dit 
elevation is at sea level, the appropriate / Kollsr 
. may be read directly from the graph shown in Fig. 1. 10 7 of the 
“ For ground-level elevations higher than sea level, the 8 Nd indica 
t appropriate correction value can be obtained by moving € a pressu 
: up the altitude scale (ordinate) to the altitude at the 4 VA for ea 
‘ ¥ ground level then moving out along this altitude value 2 each ¢ 
until the plotted curve is intersected. The abscissa re) (Temy 
of this point of intersection is the necessary correction to 20 30 40 50 6 level « 
value that must then be subtracted from the value of Pa (b) 
~~/P/T shown for the launching altitude. Thus, if cas the no 
the ground level is 1,000 ft., the value 0.15 would be Fic. 1. Graphof VP/T values, computed from U.S. Standard é otra 
: s Atmosphere Tables, for use in calculating parachute descent 
subtracted from the value for the launching altitude. times from the formula the no 
If we wish to calculate the descent time from 20,000 descent time = (302.52 © WP/T)/V/2'T/C)S altitud 


rd 
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Table for Conversion of Pressure Porometers 


ait-ft. 5600 -2000 -1000 
p-in Hg 30 320 BO 300 


© 00 2000 3000 4000 5000 6000 7000 8000 9000 
220 280 270 260 250 240 23 22 a 


20000 22000 24000 26000 28000 0000 
2 " 10 9 
30000 i 32000 i i i 36000 i 38000 i 
55 5 45 439 wv 35 


Fic. 2. The actual, or true, pressure at a given altitude is 
computed by setting one leg of a pair of dividers on zero and 
the other leg on the appropriate Koilsman number and by 
applying this correction to each altitude in question. Table 
adopted from Bellamy.” 


rule out the effect of changing atmospheric conditions 
between tests. 

For such a calculation the Kollsman number, or 
altimeter setting of the day, the temperature of the air 
at each 2,000-ft. level, the exact indicated altitude and 
temperature at time of drop, and a knowledge of the 
ground-level elevation are needed for each flight. The 
altitude of release need not be the same on successive 
flights. With these data at hand the method of calcu- 
lation is as follows: 

(a) True pressure can be obtained by entering Fig. 2 
with a pair of dividers, one leg of which is set on zero 
altitude (29.92-in. pressure) and the other leg set on 
the figure for the altimeter setting of the day.? Thus, 
if such settings were greater than 29.92, the top line 
would be used, and, if it were less than 29.92, the 
second line would be used. The setting of the dividers 
so obtained represents a pressure correction that is used 
for each 2,000-ft. level as follows: Place one leg of the 
set dividers on the 2,000-ft. level in question and swing 
the other leg to right or left depending on whether the 
Kollsman number for the setting was to right or left 
of the zero point of 29.92 in. Hg. Read the pressure 
indicated by the second leg. This represents the actual 
pressure at the given indicated altitude. By doing this 
for each 2,000-ft. level, true pressures are obtained for 
each of these levels at which temperature is known. 
(Temperature readings were taken at every 2,000-ft. 
level of indicated altitude on each flight.) 

(b) Next, the true altitude can be obtained from 
the nomogram shown in Fig. 3. To use this nomogram, 
a straightedge is placed on the right-hand portion of 
the nomogram in such a way as to connect the indicated 
altitude with its appropriate temperature. This locates 


SCALES FOR CHANGING INDICATED ALTITUDE TO TRUE 


ALTITUDE 
using Flight Temperatures 
INDICATED 
INDICATED TRUE ALTITUDE 
ALTITUDE ALTITUDE ‘TEMPERATURE Thousand 
WO E50 4° F100 
30 4 
80 $40 
20 70 
60330 7TO$35 60 
soz 58 22 
45 24 
ie -20 ° 
14 -30 
26+ 14 26 
26+ 13 4 
ate 22 -40 
eeti et 
20+10 Qe 
x H 


Fic. 3. A straightedge is placed to correct the indicated 
altitude of the right-hand column with its observed temperature. 
This locates a point on the line X. The straightedge is then 
rotated about this point so gs to connect it with the indicated 
altitude in the left-hand column. True altitude is then read at 
the point where this line intersects the ‘‘True Altitude”’ scale. 


a point on the line X. Next, the straightedge is rotated 
about that point on the line X and made to connect 
that point with the observed indicated altitude on the 
scale to the left. Then the true altitude may be read 
off the “True Altitude’ scale at the point where the 
straightedge crosses that scale. 

(c) The value for ~/P/T can now be calculated for 
each of the successive levels for which true pressures 
have been determined in step (a), whose temperatures 
were recorded during flight, and whose true altitudes 
have been determined nomographically in step (b). 
The data up to this point had best be tabulated as 
shown in Table 2. 

(d) The next step involves plotting the true alti- 
tudes vs. 1/P/T values. From this graph the mean 
~/P/T value for each 2,000-ft. stratum may then be 
read directly. A second table should now be made 
showing the mean value of ~/ P/T for each complete 
2,000-ft. stratum from ground altitude up to altitude of 
release. 

The »/P/T values of the above table may now be 
summed and used in Eq. (5). Finally, Eq. (6) will 
give the C,S value for the known weight, W, of dummy 
and parachute. 

The above procedure may be illustrated with a hypo- 
thetical case in which a 200-lb. dummy with a 24-ft. 
nylon parachute (total weight 215 Ibs.) is dropped from 
a 9,875-ft. indicated altitude, the time of descent being 
recorded as 382 sec. Altimeter setting is 30.20 and 
ground level is +300 ft. The observed temperature and 
indicated altitudes are shown in Table 2, along with the 
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TABLE 2 
Example of Calculation of CaS from Flight Data 


Indicated Tem- 
Altitude, perature, True True kdl 
Ft. °C. Pressure Altitude +/P/T 
+300 +23 296 29.86 308 0.3175 
2,000 +15 288 29.09 2,020 0.312 
4,000 +13 286 26.09 4,050 0.302 
6,000 +12 285 24.21 6,150 0.292 
8,000 +11 284 22.43 8,250 0.2805 
9,875 +9 282 20.89 10,220 0.2722 
Stratum Mean +/P/T 
300 to 2,300 0.315 
2,300—4,300 0.306 
4,300-6,300 0.296 
6,300-8,300 0.2855 
8,300-9,875 0.2163* 
1.4188 


* 1,575/2,000 of 0.2758 (mean for 8,300 to 10,300 stratum). 


From Eq. (5) 
pV2 = (302.52 X 1.4188/382)? = 1.2625 
From Eq. (6) 
CaS = 430/1.2625 = 340.6 


True Alt., Thds. ft. 


Fic. 4. Graph of true altitude vs. ~P/T for a sample 
calculation of parachute descent time. From this graph the 


mean 1/P/T value for each 2,000-ft. stratum may be read 
er ig 4 These values, together with their sum, are shown in 
able 2. 


entire calculation. The top stratum will not, in general, 
be a full 2,000-ft. stratum; therefore, since the coeffi- 
cient 302.52 in Eq. (5) is based on 2,000-ft. strata, it is 
clear that the average +/ P/T value for the top stratum 
must not be added directly along with the \/P/T values 
for the full 2,000-ft. strata but must be weighted by a 
fractional factor given by top-stratum-thickness-in-feet 
per 2,000 ft. An example of this weighting process is 


included in the sample calculations of Table 2. The 
graph described in step (d) as plotted for this caleu- 
lation is shown in Fig. 4. 


LANDING VELOCITY 


The velocity with which men of different weight and 
wearing different sizes and kinds of parachutes hit the 
ground is sometimes useful information. The effect of 
carrying heavy equipment, such as paratroops are 
sometimes required to do, and the effect of high ground 
elevations can easily be calculated. 

The velocity at any altitude can be calculated from 
the data contained above by using the mean C,S values 
determined for each type of parachute. Thus, for 28-ft. 
nylon parachutes, 


= 483 
But 
pV? = 2W/C.S 


V = V2W/C.Sp = V2W/483p (11) 


The quantity p|”, as expressed above, is seen to be a 
constant during any given descent (fixed W and para- 
chute type). Therefore, p:Vi? = p:V2?, so 


V2 = V (12) 


where 
V, = velocity at first altitude 
V2 = velocity at second altitude 
pi/p2 = ratio between the densities of the two 


altitudes 


The validity of the relationship expressed above can 
be put to test using the observed descent times listed 
in Table 1. In Fig. 5 the observations on descent times 
of 28-ft. nylon parachutes have been averaged for each 
of four approximate altitudes and plotted as launching 
altitude versus total time of descent to 2,800 ft. (ground 
level in these tests). Then a theoretical curve has been 
constructed using Eq. (12), where V2 = velocity at 
2,800 ft., pp = density at 2,800 ft., and the resulting 
curve is plotted to the same axes. The V2 used in this 
calculation was derived from Eq. (11), using standard 
p for 2,800 ft. 

It is evident that the observed mean descent times 
from the four different altitudes fall very close to this 
theoretical curve. Since the same relative correlation 
was found when the 24-ft. nylon and 28-ft. silk para- 
chute descent time data were plotted against such a 
theoretical curve, it would appear to be safe to accept 
the validity of Eq. (12). 

If we wish to determine the sea-level velocity of men 
weighing 225 lbs. wearing a 28-ft. nylon parachute, we 
may refer to Eq. (11) above and use any convenient 
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Density Ait, Thds. ft. 


l 1 l l 


200 400 600 800 1000 1200 1400 


Time in Seconds to 2800 ft. Ground Level 


Fic. 5. Descent times of 28-ft. nylon parachutes to 2,800-ft. 
ground level. The curve is derived from the formula 


Va = Vo.s00/V pa/p2,s00 


where Va = velocity at launching altitude, V2,s00 = velocity at 
2,800 ft., pa = air density at launching altitude and pe2,s00 = air 
density at 2,800 ft. The four points represent means of four 
groups of observed descent times. 


altitude. If we choose 10,000 ft., p: becomes 0.001756 
and V; becomes 23.11. Then, using Eq. (12), 


23.11 X 0.859 = 19.85 ft. per sec. at sea level 


where 0.859 is the square root of the density ratio for 
10,000 ft. The comparable sea-level velocities for 24-ft. 
nylon and 28-ft. silk parachutes are 23.78 and 18.96, 
respectively. 

The precise effect of weight change on the rate of 
descent of parachutes cannot be shown by our data, 
since we used only one weight dummy. In general, we 
know that the heavier the weight the faster the rate of 
descent and, hence, the greater the landing impact; 
but one fundamental question remains unanswered: 
Does the weight of the incorporated air mass combine 
with the effect of dummy weight to determine the 


TABLE 3 
Calculated Effect of Varying Weights of Men on the Sea-Level 
Velocity of a Standard 24-Ft. Nylon Parachute 


Velocity, 
Weight, Sea Level, 
Lbs. Ft. per Sec. 
150 19.38 
175 20.93 
200 22.37 
225 23.78 
250 25.02 
275 26.27 
300 27.40 


actual descent rates? If it does not, then the effect of 
the weight of the man on descent rate is easily calcu- 
lated as follows: 

For 24-ft. nylon parachutes with man-shaped dum- 
mies, the C,.S value is 336. Then, 


pV? = 2W/C.S = 2W/336 


Next, applying the above in Eq. (11), using the appro- 
priate value of p in slugs per cu.ft., the velocity may be 
calculated. Table 3 shows the results of such a calcu- 
lation for sea-level velocities. If the weight of the air 
mass does enter into the problem, it would have a damp- 
ing effect on the change in velocity, since the incorpo- 
rated air mass should remain relatively constant with 
changes in weight of men. Hence, the percentage 
change in weight of the whole system would be less and 
the total spread in velocity shown in Table 3 would be 
too great. Thus, for the larger weights in Table 3, the 
velocities may involve small errors, but it is to be 
noted that any errors are on the safe side, since, at 
worst, the landing velocities in Table 3 are overesti- 
mated. 

In this treatise no account is taken of the effect of 
the so-called “‘ground reflection”’ of air on the landing 
velocity of the parachute. This phenomenon, said by 
Stasevich? to give a retardation of 10 to 20 per cent, isa 
variable factor and has not, to our knowledge, been 
subjected to scientific measurement. This effect is 
presumed to be due to two factors: (1) rising air cur- 
rents near the ground due to ground heat and (2) the 
striking of the ground of the air mass being pushed 
ahead of the canopy. 

Practical experience in both the United States and 
British airborne services and air forces has shown that 
an increase in ground altitude of 5,000 ft. has been suffi- 
cient to increase greatly the number of landing injuries 
beyond those encountered in descents to near sea level.‘ 
It is generally assumed that this rise of injury rate is 
due to a rise in landing velocity occasioned by the 
reduced density of the air at the ground level in elevated 
To show more clearly the effect of altitude on 
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32 
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Fic. 6. The curves are derived from the formula Va = Vy + 
V pa/po, where Va = velocity at any altitude, V, = velocity at 
sea level, and pa/p) = the density ratio. 
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landing velocities, Fig. 6 represents the theoretical in- 
crease in landing velocities of 24 and 20 ft. per sec., 
respectively.® 

It may seem surprising, since the landing velocity at 
5,000 ft. in a 24-ft. nylon parachute is less than 2 ft. 
per sec. faster, that the injury rate should be appreci- 
ably higher. However, it must be recalled that the 
kinetic energy of a falling body increases as the square 
of the velocity. This relationship is given by the 
formula 


E, = '/.MV? 


Thus, calculation shows that the kinetic energy to be 
dissipated upon striking the ground at 5,000 ft. is some 
17 per cent greater than at sea level for the same weight 
of man and same type of parachute. On the other 
hand, it may be calculated from the same formula that 
the kinetic energy associated with the 28-ft. parachutes 
is considerably less than for the 24-ft. parachutes at 


1947 


the same ground level. In the case of the 28-ft. nylon, 
for a 200-Ib. man, it is some 43 per cent less, and in the 
case of the 28-ft. silk, it is 57 per cent less at sea level, 
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Letter to 


Dear Sir: 

Mr. Krzywoblocki’s Letter to the Editor (JOURNAL OF THE 
AERONAUTICAL SCIENCES, pp. 23-24, January, 1947) is so sig- 
nificant and representative of the opinion of a certain group of 
students of the problem of jet and rocket propulsion that it 
seems desirable to discuss further the definition of propulsive 
efficiency. 

The formula for propulsive efficiency was established correctly 
many decades ago. It reads, using Mr. Krzywoblocki’s nota- 


tion: 
2[ (Wa + Wy)x; WaVIV a; 
[(We + — WaV?] 


Let us remember a few basic laws and principles of elementary 


physics: 
(1) Velocity and force are vectors; the same is true for momen- 


tum, while kinetic energy is a scalar. Therefore, the direction of 
velocity vector is of importance in computing momentum, while 
it is immaterial in computing kinetic energy. 

(2) The principle of relativity, a few centuries old, still holds 
true in classic mechanics. 

(3) Principle of conservation of energy and principle of con- 
servation of momentum invariably hold true. 

(4) Every problem of mechanics must be placed in a certain 
given system of reference before it is theoretically treated. Of 
course, there is no difference as to the system chosen—i.e., one 
system of reference is as good as another. But once chosen, it 
must be respected until the whole computation of either energy or 
momentum balance is accomplished. Leaping from one system 
to another during the computation is inadmissible. 

(5) In the very nature of the notion of efficiency lies the fact 
that its 100 per cent value is a limit that is either asymptotic in 
function of at least one of the physical parameters involved or 
occurs for zero value of at least one other physical parameter in- 
volved, which is essentially positive (i.e., absolute temperature, 
This limit is, of course, not actually attainable, 


mass, etc.). 
This means that practically 


but it must be logically admissible. 


attainable values of efficiency may be imagined as close to 100 


the Editor 


per cent as the technical facilities now available permit, and 
there is no limit in improving these facilities so as to approach 
more closely the 100 per cent value in the future. Any formula 
for efficiency giving a theoretical maximum value of efficiency 
less than 100 per cent is erroneous. To illustrate the above, 
let us take the theoretical efficiency of the Otto cycle, which has 
an asymptotic limit of 100 per cent for the volumetric compression 


ratio infinitely large. Another example is the efficiency of Carnot 
cycle 
_ — 
n 


which has either an asymptotic limit of 100 per cent for the abso- 
lute temperature 7; infinitely large or a limit of 100 per cent for 
the absolute temperature 72 equaling zero; of course, the abso- 
lute temperature cannot be negative. 

Mr. Krzywoblocki’s assumption (2) is erroneous. Taking the 
propelled unit as a system of reference, air enters the engine at 
the rate of Wz lbs. per sec. and with the speed V; the input 
momentum is, therefore, (W./g) V and the input kinetic energy is 
(W,/2g) V?; this is an actual fact. Ascribing this latter energy 
to the thermal efficiency would only displace its asymptotic limit 
from the 100 per cent value, which is wrong. The energy to be 
produced by the engine is, therefore, only 


[((Wa + W,)/2g] — (Wa/2g) V? 


Now, Mr. Krzywoblocki changes to another system of refer- 
ence. Energy L; given by his Eq. (2) was computed relative to 
the propelled unit. Now he adds the kinetic energy Lz, given by 
his Eq. (3), relative to the earth. This energy is, of course, nil 
in the former system of reference. A further step along the path 
chosen by Mr. Krzywoblocki may be made. The earth as a sys- 
tem of reference is no better than any other system. Let us take 
the sun as a system of reference. The velocity V in Mr. Krzywo- 
blocki’s Eq. (3) for Lz will become about 90,000 ft. per sec. The 
‘“‘most favorable conditions’ in which the maximum energy can 
be imparted to the jet unit become still more favorable. Any 
other fixed star or even the Milky Way may also be chosen as 4 
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system of reference. The only condition is to keep the direction 
of fight perpendicularly to the corresponding gravitational field. 
But this is not all. The internal energy of the fuel (not to be 
confused with the intramolecular—i.e., chemical energy) is, in 
fact, molecular kinetic energy. Why not add this energy also? 
True, molecular velocities have many different directions, but this 
makes no difference, since kinetic energy is not a vector. Men- 
tion should be made of the intra-atomic energy also. 

Mr. Krzywoblocki’s statement that ‘“‘undoubtedly, the best 
speed of an airplane occurs when the optimum value of the useful 
power may be achieved”’ is misleading. It follows from his fur- 
ther reasoning that he speaks about the unit power, as obtained 
per 1 lb. of air-involved. One B.t.u. in hydrogen weighs three 
times less than in gasoline. Does this mean that hydrogen is 
three times better than gasoline as fuel for combustion engines? 
Of course, the increase in air output of a jet engine will increase 
its bulk. But we are discussing only the problem of efficiency 
of the engine, other aspects of the propelled unit being, there- 
fore, logically immaterial here. But even if they were involved, 
the best total useful power—and, consequently, the best speed of 
an airplane—would certainly not occur at the maximum of the 
unit power but somewhere between this maximum and the maxi- 
mum efficiency. ; 

The same kind of misunderstanding appears in Mr. Krzy- 
woblocki’s criticism expressed as follows: ‘‘One may ask what 
this ‘propulsive’ efficiency represents if it is not equal to zero and 
both the propulsive force and propulsive power are equal to 
zero.” This applies to the “‘airscrew” case, W;/W, ~ 0, thus 
ny = 2V/(v; + V), according to Eq. (1), where v; = V. The 
point is that the statement that propulsive force and power are 
zero here is untrue. The air output is in this case infinitely large 
and the product (© X 0) is finite. This is, of course, an abstrac- 
tion. But similar abstractions occur in theoretical physics so 
often that every physicist and even every college student must be 
familiar with them. The abstractive limit (© XX 0) cannot be, 
of course, actually achieved. But we can approach it as close as 
the actual technical facilities now available permit. A very large 
output together with an extremely small increase in velocity may 
be actually applied, thus giving an efficiency close to the theoreti- 
cal limit. The usual airscrew is certainly close to this limit. A 
100 per cent efficiency is here an asymptotic limit in function of 
air output W,. Assuming the propulsive force is given as (Th), 
we get: 


(2) 


V + [g(Th)/WaVj \ 
nm = 2V/{2V + [g(Th)/Wa]} 


Let us, finally, consider the case of a pure rocket, where W, = 0 
in Eq. (1), i.e.: 


Wy; 
propulsive force (Th) = 
propulsive power P = V(Th) = ——— (3) 
2V 2V 
efficiency np = = 


where w = v,; — V is the velocity of burnt gases relative to the 


earth. Consider the earth as a system of reference. During 
the energy release (i.e., by combustion) the velocity of the fuel 
relative to the earth changed from (—V) tow. Therefore, the 
energy released, if its total is represented as kinetic energy (ther- 
mal efficiency 100 per cent), will correspond to the numerical 
value of the vecto~ [w — (— V)], i-e.: 
ay, 2 
Pus, =o! (w — (— = 
2g 2g 
The fuel traveling with the rocket had, before combustion and 
relative to the earth, kinetic energy (W;/2g)V?. The burnt gases 
leave the engine with the kinetic energy (W;/2g)w?, relative to 
the earth. The energy lost, which is gone uselessly to the sur- 
roundings, is therefore 


Pioss = (W;/2g) (w* — V*) (4) 


and this energy has been borrowed from the total energy released. 
The fact that V and w have opposite directions is immaterial, 
since kinetic energy is a scalar. Therefore, the useful energy 
must be 


Prot.—Pioss = (W;/g) Vo; 


which is in accordance with the propulsive power P found directly 
as a product of velocity V and force (Th) acting with this velocity 
(always relatively to the earth). 

The energy lost Pioss as expressed by Eq. (4), is an essentially 
positive physical quantity. It would be illogical to suppose that 
this loss, dispersed in the surroundings, may be negative—i.e., 
that kinetic energy may be derived from surroundings that are at 
rest. In other words, it is impossible to transfer any part of the 
initial kinetic energy of fuel W;V?/2g to the rocket, since both 
move with the same speed V—i.e., one is at rest relative to an- 


other. The theoretical upper limit for efficiency is, therefore, 
given by Pioss > 0, i.e.: 
(5) 


Impulse drawn from the reaction of surroundings does not exist 
here. The picture is different when it does—e.g., in a gun, where 
the impulse balanced by the reaction of the earth is imparted and 
the theoretical limit for velocity becomes V < vj. The case of 
pure jet (‘‘airscrew’’) is somewhat analogous. Here, the impulse 
of surroundings (which belong to the earth system and are at rest 
relative to it) exists and may be explained as being in balance 
with the reaction of the surrounding air; this reaction is due to 
the inertia of air mass. In the general case, where W, ~ 0 and 
W; # 0, the impulse exists and involves in the energy balance a 
certain work. 

The new ‘‘theories’’ as criticized above may cause so much 
harm in the era where the problem of dynamic propulsion is so 
intensively studied and so much time and money has been and still 
is being spent for experiment that it seems desirable to refute 
them at once. . 

B. SzczZENIOWSKI 
Associate Professor 
Ecole Polytechnique 

University of Montreal 
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Services 
of the Libraries of the 


Institute of the Aeronautical Sciences 


The services of the Libraries are available to all mem- 
bers of the Institute, to Corporate Members, to advertisers 
in the AgronauticaL ENGINEERING Review and 
NAUTICAL ENGINEERING CaTALoa, and, under usual library 
limitations, to the public. Four specialized services are 
available. 


The Paul Kollsman Lending Library 


This lending library service makes available, without 
charge, the latest and more important aeronautical books. 

Members may request the loan of any aeronautical or 
technical book they wish to borrow. Through an exchange 
agreement with the Engineering Societies Library, any 
book on general engineering may be borrowed from its great 
collection of over 160,000 volumes. 

A photostating service is available at usual library rates. 

Applications for membership in the library and further 
information will be sent on request. 


The W. A. M. Burden Reference Library 


This reference library contains over 12,000 aeronautical 
books, magazines, pamphlets, and reports gathered from 
world-wide sources and is one of the most complete aero- 
nautical libraries in the world. Material from this library 
is not available for loan but may be used for reference pur- 


poses. 


The Pacific Aeronautical Library 


6715 Hollywood Boulevard 
Los Angeles 28, California 


Established in cooperation with the aircraft companies 
the library serves. The leading aircraft companies in or 
near Los Angeles participate in its support and operation. 

This service library for aeronautical research is available 
to the public for reading privileges. Source material in- 
cludes aerodynamic and structural research reports, as well 


INSTITUTE OF THE AERONAUTICAL SCIENCES, INC. 
2 EAST 64TH STREET, NEW YORK 21, N.Y. 


as books on drafting, production methods, history, and al- | 


lied sciences. It furnishes books, periodicals, and pamphlet | 
material to the participating aircraft companies to supple | 


ment their engineering libraries. 


Technical Information Service 


This service has experienced personnel under the super- 
vision of trained aeronautical engineers to compile any in- 
formation desired. The services range from listing special- 
ized reference books to the preparation of exhaustive 
bibliographies, digesting of reports, and general surveys of 
any acronautical subject. Some of the available services 
are: 


Bibliographies on any aeronautical subject. 

Reports on any aeronautical subject. 

Digests of aeronautical books, papers, periodicals, and refer- 
ences. 

Translations. 

Engineering investigations of special aeronautical subjects. | 

Biographies of individuals engaged in aeronautics. 

Photostats of any aeronautical or general enginecring 
material. | 

Microfilms made on special order. 

Photographs made from the Institute's photographic coilec- 
tion. 

Drawings and tracings made. 


In addition to the services'mentioned any commission 
which comes within the scope of the Service will be ac 
cepted. Special arrangements may be made for work re- 
quiring several weeks or months. 

Translators are available for accurate transcriptions of all 
foreign language data. Translations are carefully edited 
by trained engineers. 

Reproductions of any material in the collections of 
the Institute may be ordered at standard photostat rates. 
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